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I.   NOMENCLATURE 

X   -  coordinate  in  the  fore  and  aft  direction,  X^  direction. 

Y   -  coordinate  in  the  transverse  direction,  X2  direction. 

Z   -  coordinate  in  the  verticle  direction,  X3  direction. 

M^  -  mass  of  the  left  front  wheel  and  tire. 

M2  -  mass  of  the  right  front  wheel  and  tire. 

M3  -  mass  of  the  inner  left  rear  wheel  and  tire. 

M4  -  mass  of  the  inner  right  rear  wheel  and  tire. 

M5  -  mass  of  the  main  chassis. 

Mg  -  mass  of  the  front  axle  assembly. 

M7  -  mass  of  the  outer  left  rear  wheel  and  tire. 

Mg  -  mass  of  the  outer  right  rear  wheel  and  tire. 

J^j  -  rotational  moment  of  inertia  of  M^  about  the  j  axis,   i  ■  1,  2,  ...  8 

j  =  X,  t,    Z. 
Kjj  -  spring  rate  of  Mj_   in  the  Xj  direction,   i  =  1,  2,  3,  4,  7,  8; 

j  =  1,  2,  3. 
Cjj  -  damping  coefficient  of  Kj_   in  the  Xj  direction,   i  =  1,  2,  3,  4,  7,  8; 

i  -  1.  2.  3- 
y-jj  -  ground  profile  seen  by  the  tire  to  ground  contact  point  of  M^  in  the 

Xj  direction,  i  =  1,  2,  3,  4,  7,  8;  j  =  1,  2,  3. 
y^j  -  the  first  derivative  of  yjj  with  respect  to  time. 
EI  -  the  product  of  the  modulus  of  elasticity  in  tension  and  the  cross 

sectional  moment  of  inertia  of  the  rear  axle. 
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JG  -  the  product  of  the  cross  sectional  polar  moment  of  Inertia  and  the 

modulus  of  elasticity  in  shear  of  the  rear  axle. 
x±     -  motion  coordinates,   i  =  1,  2,  3,  ...  29. 
x-l     -   the  first  derivative  of  x^  with  respect  to  time. 
8^  -  the  second  derivative  of  x^  with  respect  to  time. 
AXi  -  x± 
VXi  -  Xl 


II.   INTRODUCTION 

The  design  of  large  off  road  vehicles,  such  as  farm  tractors,  is  a 
process  which  makes  use  of  many  different  methods  to  determine  vehicle 
parameters  and  predict  vehicle  performance.   One  such  technique  is  the  use 
of  mathematical  modeling  and  the  digital  computer  to  simulate  the  motion  of 
a  vehicle  in  question. 

Computer  simulations  allow  the  design  engineer  to  vary  certain  param- 
eters, and  observe  the  resultant  motion.  This  procedure  avoids  the  cost  of 
building  a  prototype,  and  modifying  that  prototype  until  the  desired  per- 
formance is  achieved.   It  also  alleviates  the  risk  of  persons  being  injured 
and  equipment  being  damaged  while  observing  a  system  that  may  fail. 

With  a  predetermined  set  of  parameters,  mathematical  modeling  may  be 
utilized  to  observe  the  performance  of  a  vehicle  subject  to  varying  surface 
conditions.   These  surface  conditions  may  be  due  to  a  certain  geometrical 
configuration,  such  as  a  side  slope  or  bump,  or  different  soil  types,  such 
as  sand  or  hard  packed  soil.  Vibrations  induced  by  tire  tread  patterns  can 
also  be  observed. 

Many  mathematical  models  of  farm  tractors  have  been  developed.  Most  of 
these  were  derived  to  describe  the  motion  of  single  rear  wheeled  vehicles. 
The  rear  axles  have  usually  been  treated  as  rigid  members,  which  greatly 
reduces  the  degrees  of  freedom  of  the  system  and  hence  the  number  of  equa- 
tions of  motion. 


In  this  thesis,  a  dual  rear  wheeled  farm  tractor  is  modeled  as  eight 
lumped  masses.  The  rear  axles  are  considered  to  be  elastic,  allowing  the 
vehicle  twenty-nine  degrees  of  freedom.  The  equations  of  motion  for  the 
system  are  derived  by  two  different  methods,  allowing  one  set  of  equations 
to  be  used  as  a  check  for  the  other. 

The  first  of  the  sets  is  generated  using  Newton's  laws  and  elementary 
beam  theory.   The  second  set  uses  the  energy  method,  taking  advantage  of 
the  Lagrange  equation  for  nonconservative  systems.   In  both  cases,  the 
equations  are  reduced  by  the  Gauss-Jordon  method  to  yield  twenty-nine  second 
order  differential  equations  of  motion.   These  differential  equations  are 
numerically  integrated  using  the  fourth  order  Runge-Kutta-Gill  method  to 
give  displacements  as  functions  of  time. 


III.   REVIEW  OF  LITERATURE 

Among  the  first  work  in  tractor  modeling  was  that  of  McKibben  [7],  done 
in  1927.   In  a  series  of  articles,  he  investigated  the  kinematics  and  dynam- 
ics of  single  rear  wheeled  farm  tractors.  McKibben  described  the  forces 
acting  on  the  vehicle  due  to  various  causes,  such  as  soil  reactions,  drawbar 
pull,  and  the  gyroscopic  action  of  the  engine  flywheel. 

These  articles  considered  the  general  motion  of  a  tractor  as  well  as  a 
discussion  on  backward  overturning  while  moving  up  a  slope.  Unfortunately, 
McKibben 's  work  was  done  with  vehicles  riding  on  steel  wheels  rather  than 
pnumatic  tires.  He  treated  the  tractor  as  a  single  rigid  mass,  and  avoided 
the  use  of  second  order  differential  equations  by  considering  only  constant 
velocities. 

Since  McKibben's  time,  tractors  have  become  somewhat  more  sophisticated, 
and  with  the  advent  of  high  speed  computers,  larger  more  complex  equations 
can  be  handled.  The  more  recent  work  in  tractor  modeling  can  be  divided 
into  two  general  categories.   The  first  type  is  concerned  with  gross  motion 
of  a  vehicle,  such  as  backward  overturning,  or  sideward  overturns  on  side 
slopes.  The  second  type  considers  small  amplitude  vibrations  for  studies 
in  rider  comfort  or  fatigue  of  tractor  components. 

A.   Studies  in  Gross  Motion 

In  1967,  Goering  and  Buchele  [2]  developed  a  model  to  predict  backward 
overturning  of  an  unsprung  vehicle.   They  restricted  the  motion  of  the  vehicle 


to  a  plane  normal  to  the  ground,  and  considered  the  rear  axles  to  be  rigid. 
The  result  of  their  work  was  a  set  of  nineteen  equations  that  were  used  to 
simulate  the  effect  of  a  sudden  clutch  engagement  with  no  drawbar  load. 

Hudson  et  al.  [5]  formulated  a  two  dimensional  stability  model,  in  1973, 
to  determine  the  behavior  of  single  rear  wheeled  tractors  on  slopes.  They 
developed  two  differential  equations  that  could  be  easily  programmed  into  an 
analog  computer.  These  equations  determined  the  angle  of  pitch  of  a  tractor, 
given  the  drawbar  load. 

A  model  to  simulate  both  backward  and  sideward  overtumings  was  devel- 
oped by  Davis  et  al.  [l]  in  1974.   Their  model  treated  a  single  rear  wheeled 
tractor  as  five  masses  with  ten  degrees  of  freedom.  The  rear  axles  were 
assumed  to  be  rigid,  so  that  the  rear  wheels  were  constrained  to  have  the 
same  motion  as  the  main  chassis,  except  for  their  ability  to  rotate  about 
the  longitudinal  axis  of  the  axle.   This  model  was  capable  of  simulating 
tractors  with  either  wide  front  ends,  or  tricycle  type  front  ends. 

B.   Studies  in  Small  Amplitude  Vibrations 

A  model  to  predict  small  amplitude  vibrations  of  a  single  rear  wheeled 
farm  tractor  was  developed  by  Raney  et  al.  [ll]  in  1961.  For  their  model, 
the  frame  and  rear  axle  were  assumed  to  be  a  single  rigid  body.  This  model 
had  only  three  degrees  of  freedom,  verticle  displacement,  pitch  and  roll.  A 
set  of  three  second  order  differential  equations  was  derived  and  programmed 
into  an  analog  computer.  The  authors  also  considered  the  effect  of  adding 
a  vibration  absorber  to  the  system. 

In  1964,  Huang  et  al.  [4]  considered  the  effect  of  elastic  rims  and 
elastic  suspension  on  the  motion  of  a  single  rear  wheeled  vehicle.  They 


derived  a  model  having  four  degrees  of  freedom,  with  which  they  could  compare 
the  effect  of  elastic  rims  to  rigid  rims. 

Pershing  and  Yoerger  [9]  formulated  a  model  for  a  single  rear  wheeled 
tractor  with  a  side  mounted  implement.  This  work,  done  in  1968,  was  used  to 
study  the  motion  of  a  tractor  while  mowing  side  slopes.  Although  the  equa- 
tions were  derived  for  small  amplitude  vibrations,  they  were  used  to  determine 
the  size  of  a  bump  that  would  cause  a  sideward  overturn  of  the  vehicle. 
Pershing  and  Yoerger  considered  the  tractor  as  five  masses  with  nine  degrees 
of  freedom,  and  treated  the  rear  axles  as  rigid  members.  The  Lagrange  equa- 
tion was  used  to  formulate  the  nine  differential  equations  of  motion,  which 
were  then  integrated  numerically  using  a  digital  computer. 

Mather  [6]  developed  a  model  in  1970  to  determine  the  stress  variations 
in  the  rear  axle  of  a  dual  rear  wheeled  tractor.  He  modeled  the  tractor  as 
eight  lumped  masses  with  nineteen  degrees  of  freedom,  and  treated  the  rear 
axles  as  elastic  cantilever  beams.   In  his  formulation,  the  Lagrange  equation 
was  used  to  derive  the  equations  of  motion. 

The  effect  of  the  front  wheels  being  allowed  to  rotate  about  the  front 
axle  in  Pershing  and  Yoerger' s  model  was  neglected  in  a  model  developed  by 
Wolken  and  Yoerger  [l3]  in  1974.  Their  work  considered  the  dynamic  response 
of  a  single  rear  wheeled  tractor  subject  to  random  inputs.  They  modeled  the 
tractor  as  five  masses  with  seven  degrees  of  freedom.  The  effects  of  forward 
velocity,  tire  spring  rate,  and  frame  moments  of  inertia  on  the  motion  of 
the  tractor  were  all  investigated. 

During  the  same  year,  Smith  and  Yoerger  [l2]  developed  a  similar  model 
to  predict  variations  in  the  forward  motion  of  a  single  rear  wheeled  tractor 
due  to  a  changing  drawbar  load.  They  also  modeled  the  tractor  as  five  masses 


with  seven  degrees  of  freedom,  but  motion  was  restricted  to  a  plane  normal  to 
the  ground. 


IV.   THE  SYSTEM  BEING  MODELED 

A.  Assumptions 

1.  A  dual  rear  wheeled  tractor  is  treated  as  eight  discrete  masses; 
four  rear  wheels,  two  front  wheels,  the  front  axle  assembly,  and 
the  main  chassis.  Each  mass  is  assumed  to  be  concentrated  at  the 
center  of  gravity  of  its  respective  body. 

2.  Small  angular  displacements  are  assumed,  to  keep  the  model  linear. 

3.  The  sections  of  the  rear  axle  between  the  frame  and  inner  wheels, 
and  between  inner  and  outer  wheels  are  assumed  to  act  as  elastic 
cantilever  beams,  with  negligable  mass  and  internal  damping. 
Longitudinal  vibrations  of  these  sections  are  neglected. 

4.  Each  tire  is  assumed  to  act  as  a  set  of  three  mutually  perpendicular 
linear  springs  with  viscous  damping. 

5.  The  reaction  between  a  tire  and  the  ground  is  assumed  to  act  at  a 
single  point. 

6.  The  tractor  oscillates  about  its  equilibrium  position. 

B.  The  Coordinate  System  and  Degrees  of  Freedom 

A  right  hand  coordinate  system  is  used  as  a  reference,  with  its  origin 

in  the  plane  of  the  ground  and  directly  below  the  center  of  the  rear  axle. 

The  positive  X  axis  is  in  the  forward  direction  of  travel  of  the  tractor, 

the  positive  Y  axis  is  to  the  right  as  viewed  by  an  observer  riding  on  the 
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tractor,  and  the  positive  Z  axis  is  toward  the  center  of  the  Earth.  Positive 
rotations  are  taken  according  to  the  right  hand  rule. 

The  main  chassis  is  free  to  translate  in  any  of  the  three  coordinate 
directions,  and  to  rotate  about  any  of  the  three  coordinate  axes.  Hence, 
it  has  six  degrees  of  freedom.  Because  of  the  third  assumption,  the  rear 
wheels  are  constrained  to  move  in  the  Y  direction  with  the  main  chassis. 
They  are,  however,  allowed  to  translate  in  the  X  and  Z  directions,  and  to 
rotate  about  any  of  the  three  axes,  giving  each  rear  wheel  five  degrees  of 
freedom. 

The  front  axle  assembly  is  constrained  such  that  it  has  no  motion  with 
respect  to  the  main  chassis,  except  for  rotation  about  its  pin  connection 
with  the  chassis.  The  front  wheels  are  constrained  such  that  they  have  no 
motion  with  respect  to  the  front  axle  assembly,  except  for  rotation  about 
the  longitudinal  axis  of  the  axle.  Each  of  these  three  bodies  have  only 
one  degree  of  freedom. 

As  a  unit,  the  tractor  then  has  twenty-nine  degrees  of  freedom. 
Figure  1  shows  these  variables  along  with  the  reference  axes. 
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V.   EQUATIONS  OF  MOTION 

A.  Derivation  Using  Newton's  Laws 

The  equations  of  motion  are  first  derived  using  Newton's  laws  and 
elementary  beam  theory.  To  do  this,  the  forces  acting  on  the  system  must 
be  considered.   These  forces  are  at  the  tire  to  ground  contact  points,  and 
the  points  at  which  wheels  join  axles  or  axles  join  the  frame.  There  are 
seventeen  such  points,  six  tire  to  ground  and  eleven  member  to  member. 

At  each  point,  the  force  present  is  broken  into  components  in  the  X, 
Y  and  Z  directions,  and  the  moment  is  broken  into  components  about  each 
coordinate  axis.  There  are  three  force  components  at  each  point  accounting 
for  fifty-one  unknown  forces.  The  same  is  not  true  of  the  moments  at  each 
point,  since  some  members  are  free  to  rotate  with  respect  to  other  members. 

The  front  wheels  are  free  to  rotate  about  the  front  axle;  therefore, 
no  moments  are  present  at  these  points  in  the  xg  and  xg  directions.  The 
front  axle  assembly  is  allowed  to  rotate  about  its  pin  connection  with  the 
frame,  in  the  X7  direction,  eliminating  the  possibility  of  a  moment  at  this 
point  in  that  direction.  At  the  tire  to  ground  contact  points,  it  is  assumed 
that  there  is  no  moment  present  in  any  direction.   There  are  three  moment 
components  at  every  other  point,  for  a  total  of  thirty  unknown  moments. 
Figures  2,  3,  4  and  5  show  the  forces  and  moments  present  at  each  point,  and 
their  assumed  directions. 
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The  system  may  be  forced  by  a  drawbar  load,  or  by  varying  the  surface 
over  which  the  tractor  moves.   The  first  type  of  forcing  function  is  broken 
into  components  in  the  X,  Y  and  Z  directions,  and  input  directly  into  the 
dynamic  equations.  The  second  type  is  broken  into  the  surface  profile  that 
each  contact  point  sees  in  each  coordinate  direction.  These  horizontal, 
transverse  and  verticle  components  are  input  by  subtracting  them  from  the 
deflection  of  each  contact  point  in  the  X,  Y  and  Z  directions,  respectively. 
The  new  deflections  are  used  when  applying  Hooke's  law  to  the  tire  ground 
reaction.  When  this  method  is  used  to  force  the  model,  both  the  functions 
and  their  first  derivatives  with  respect  to  time  must  be  prescribed. 

There  are  twenty-nine  unknown  displacements,  fifty-one  unknown  forces, 
and  thirty  unknown  moments  for  a  total  of  one  hundred  ten  unknowns.   It  is 
then  necessary  to  have  one  hundred  ten  equations  to  describe  the  system. 
These  equations  come  from  four  sources,  static  equations  involving  the  forces 
on  the  rear  axles,  dynamic  equations  describing  the  motion  of  each  mass, 
elastic  equations  due  to  the  bending  of  the  rear  axles,  and  equations  involv- 
ing the  application  of  Hooke's  law  to  the  tires. 

a.   Static  Equations 

The  sections  of  the  rear  axle  between  the  frame  and  inside  wheels,  and 
between  inner  and  outer  wheels,  are  assumed  to  be  massless.  Therefore,  the 
following  equations  must  be  satisfied  for  each  of  these  sections,  in  each 
coordinate  direction 

£*i  =  0    and     ETj  =  0  (la.b 

Applying  (la)  to  the  section  of  axle  between  the  inner  and  outer  left 
rear  wheels,  in  the  X,  Y  and  Z  directions,  respectively,  yields 
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-F77  -  F34  =  0   ,   -F78  -  F35  =  0   and  -F79  -  F36  =  0       (2a,b,c 

Applying  (lb)  to  the  same  9ection  about  the  M3  end  yields 

-T77  -  T34  -  F36  •  L37  =  0   ,   -T78  -  T35  =  0   and 

"T79  "  T36  +  F34  *  L37  =  ° 


(3a,b,c 


Equations  (2a)  through  (3c)  may  be  rearranged  to  give  expressions  for 
F77>  F78>  F79>  T77>  T78  and  T79 


F77  =  -F34 

F78  "  _F35 
F79  =  -F36 


T77  "  _T34  "  F36  *  L37 

T78  =  -T35 

T79  =  _T36  +  F34  *  L37 


(4a-f 


In  a  like  manner,  the  forces  and  moments  acting  on  the  other  sections  are 
summed  to  give 


F54  =  -F 


37 


'55 


c56 


-F 


38 


-F 


39 


F57  "  "F44 


F58 

F59  " 


c45 


46 

F84  =  "F47 
F85  =  _F48 


r86 


-F 


49 


^54 


-T,7  +  F 


37  T  r39   L53 


T55  =  -T 


33 


x56 


-T 


39 


F37  '  L53 


T57  "  "T44  "  F46  *  L54 

T58  =  "T45 

T59  =  _T46  +  F44  '  L54 

T84  "  -^47  +  F4g  •  L84 

T85  "  _T48 

T86  "  "T49  "  F47  "  L84 


(5a-f 


(6a-f 


(7a-f 
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b.  Dynamic  Equations 

Similar  equations  can  be  written  by  summing  the  forces  and  moments  about 
the  center  of  mass  of  each  of  the  masses  one  through  eight.  The  equations 
that  must  now  be  satisfied  are 

ETj  ■  M±X±     and     IT±   =  J^  (8a, b 

Application  of  (8a)  and  (8b)  to  the  left  front  wheel  yields 


Fll  +  F17  ■  Ml  * 

(Xj_  +  (D53  +  D62)  •  X5  +  (D99  +  D13  +  Y)  •  Xg) 

(9 

F12  +  F18  "  Ml  • 

(*2  -  D53  •  K4  +  D51  •  X6  -  D62  •  *?) 

(10 

F13  +  F19  *  Ml  • 

(X3  -  Y  •  X4  -  D51  •  X5  -  (D99  +  D13)  •  X7) 

(11 

■F12  •  D14  +  F19 

"  D13  +  T17  =  J1X  '  K7 

(12 

Fll  -  D14  =  J1Y 

•  K8 

(13 

'F17  '  D13  +  T19 

=  J1Z  '  K6 

(14 

Likewise,  the  forces  and  moments  are  summed  about  the  center  of  mass  of 
masses  two  through  eight,  respectively,  to  give  (see  Figure  6) 

F21  +  F24  =  M2  •  <>L  +  (D53  +  D62>  '  s5  "  (D99  +  D22  -  ?)  '  *6>  d5 

F22  +  F25  =  M2  *  (X2  "  D53  '  *4  +  D51  "  H  "   D62  '  *7>  (16 

F23  +  F26  "  M2  '  (K3  "  ?  "  X4  "  D51  '  X5  +  CD99  +  D22)  '  V  (17 

-F22  •  D,4  -  F26  •  D22  +  T24  =  J2X  •  X7  (18 

F21  *  D24  =  J2Y  •  S9  (19 

F24  •  D22  +  T26  "  J2Z  '  h  <2° 


F31  +  F34  +  F37  *  M3  •  X14  <21 

F32  +  F35  +  F38  "  M3  '  (*2  "  D57  '  *4  "  *  '  *6>  (22 

F33  +  F36  +  F39  "  M3  '  X10  (23 
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F32  '  D34  +  F33  *  D31  "  F36  '  D32  +  F39  '  D33  +  T34  +  T37  =  J3X  '  x26  <24 

(25 
(26 


-F31  •  D 
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c37 


s22 


c44 


F42  +  F45  +  F48  "  M4 


x15 

(«2  -  D5?  •  «4 


X 


c42 


c41 


u44 


+  F, 


43 


11 

46 


u44 


M8 


-F41  •  D41  +  F44 


X  •  x6) 


J42 


+  F 


^9 


u43 


+  T44  +  T 


47 


J4X 


D42  "  F47  '  D43  +  T46  +  T49  "  J42  '  X23 


(27 
(28 
(29 
X27  (30 
(31 
(32 


F51  +  F54  +  F57  +  FDB1  =  M5  *  xl 
F52  +  F55  +  F58  +  m2  "  M5  '  K2 
F53  +  F56  +  F59  +  FDB3  "  M5  '  S 


-F 


52  *  u53 


'53  •  Y  -  F55  •  D57  -  F56  •  D56  -  F5g  •  D5?  +  F 


59 


D55" 
F51  •  D53  -  F53  • 

X  +  FDB,  • 


>55  -  FDB2  •  DB3  -  FDB3  •  Y  +  T 


54  +  T57  =  J5X 


?57  '  D57  +  F59 
I,  •  DB3  +  FDB3  •  DB-|_  +  T 


'53  •  D51  +  F« 


54  •  D57  +  F56  •  X  +  F 

'52  +  T55  +  T 


F51  •  Y  +  F52  •  D51  +  F54  '  D56  "  F 
X  +  FDBj^  •  Y  -  FDB2  • 


r5Y  *  X5 
57  '  D55  "  F58  " 


DB 


'55  '  X  "  F 
'l  +  T53  +  T56  +  T59  "  J5Z  ' 


:58  "  J5 
'58 
'5Z  •  *6 


(33 
(34 
(35 

(36 

(37 

(38 


(5c,  +  (D„  +  D-.)  •  k.  +  Y  •  k,) 


J53 


'63 
H  *   D51  '  X6  -  D63  •  K7> 


'51 


*,) 


(39 

(40 
(41 


F18  ' 

D64 

+  F19 
T17" 

*  D99  +  F25 
T24  =  J6X  ' 

•  D64 

h 

F26  • 

D99 

"F52  ' 

'  D63 

'F17  • 

D64 

"F24 

•  D64  +  F51 

•   D63 

- 

T52  = 

J6Y 

•  S5 

•F17  * 

D99 

+  F24 

•  D99  -  T19 

-he 

- 

T53  = 

J6Z 

•  h 

F71  + 

F77 

=  M?  • 

he 

F72  + 

F78 

=  M7  • 

«2  -  D5?  • 

*4  - 

X 

•  V 

F73  + 

F79 

=  M?  • 

X12 

■F72  * 

D74 

+  F73 

*  D71  +  F79 

•  D73 

+ 

T77  = 

J7X 

'   S28 
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(42 
(43 
(44 


(45 
(46 
(47 
(48 

F71  •  D74  +  T78  "  J7Y  '  *20  <49 

"F71  *  D71  "  F77  •  D73  +  T79  =  J7Z  "  X24  (5Q 

F81  +  F84  "  "8  '  S17  <51 

F82  +  F85  =  "8  '  »2  "  D57  •  X4  "  f  *  V  <52 

F83  +  F86  =  "s  •  *13  (53 

"F82  '  D84  +  F83  *  D81  "  F86  '  D82  +  T84  =  J8X  '  K29                 (54 

F81  •  D84  +  T85  =  J8Y  '  X21  <55 

"F81  •  D81  +  F84  '  D82  +  T86  =  J8Z  '  X25  <56 

c.  Elastic  Equations 

Using  the  rules  of  elementary  beam  theory,  and  the  principle  of  super- 
position, a  set  of  equations  is  written  to  describe  the  deflection  and  slope 
of  the  rear  axles  due  to  the  forces  and  moments  applied  to  them.   Consider 
the  deflection  of  the  left  rear  axle  in  the  vertical  direction.  Figure  7 
shows  this  case.  The  equations  for  the  deflections  and  slopes  of  the  inner 
and  outer  wheels  in  this  plane  are 


23 


24 


-F39  •  (L533/3EI)  +  T37  •  (L532/2EI)  =  -x3  +  (D5g  +  L53)  •  x4  - 

X  •  X;  +  X, g  +  D33  '  Xjg  (57 

F39  "  (L532/2EI)  -  T37  •  (L53/EI)  =  -x4  +  x26  (58 


and 

-F?9  •  (L3?3/3EI)  +  T?7  •  (L372/2EI)  =  -x1Q  +  x12  +  (D32  +  L3?)  • 

x26  +  D73  •  x28  (59 

F79  '  (L372/2ED  -  T77  •  (L37/EI)  -  -x26  +  x28  (60 

Similar  sets  of  equations  are  written  for  the  deflection  and  slope  of  the 
left  axle  in  the  horizontal  plane,  and  for  the  right  axle  in  the  verticle  and 
horizontal  planes. 

-F37  •  (L533/3EI)  -  T39  •  (L532/2EI)  =  -Xl  -  D57  •  x5  -  (Djg  +  L53)  • 

x6  +  x14  "  D33  "  x22  (61 

-F37  •  (L532/2EI)  -  T39  •  (L53/EI)  -  -Xg  +  x22  (62 


-F77  •  (L373/3EI)  -  T79  •  (L372/2EI)  -  -xu  +  x16  -  (D32  +  L37)  • 

x22  -  D;3  •  x24  (63 

-F77  •  (L372/2EI)  -  T79  •  (L37/EI)  =  -x22  +  x24  (64 


"F46  -  (L543/3EI)  -  T44  •  (L542/2EI)  =  -x3  -  (D55  +  L54)  •  x4  - 


(65 
46  -  U-J4-/«U  "  ^44  "  ^54' Bi'  "  "x4  +  x27 


-F46  •  tts42/2EI)  -  T44  •  (L,4/EI)  =  -x,  +  x„  (66 
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"F86  *  (L843/3EI>  "  T84  *  a842/2EI)  =  -xu  +  x^  -  (D^  +  I_84)  • 

x27  —   82  *  ^29  (o7 

"F86  •  <L842/2EI)  "  T84  •  <L84/EI)  "  "x27  +  x29  <68 

"F44  '  (L543/3EI)  +  T46  •  (L542/2EI)  =  -3^  -  D5?  •  x5  +  (D55  +  L54)  • 

X6  +  X15  +  D42  '  X23  (69 

F44  *  (L542/2EI)  -  T46  •  (L54/EI)  -  -xg  +  x23  (70 

-F84  •  (L843/3EI)  +  T8fi  •  (L842/2EI)  =  -Sy  +  x„  +  (D43  +  Lg4)  • 

F84  *  (L842/2ED  -  T8fi  •  (Lg4/EI)  =  -x23  +  x25  (72 

Elementary  beam  theory  Is  also  used  to  determine  the  amount  of  twist  in 
the  rear  axles.  The  angle  of  twist  in  a  round  bar  of  length  L  subject  to  a 
moment  T  is  given  by 

8  =  TL/JG  (73 

Applying  (73)  to  the  sections  of  the  rear  axle  yields 

T38  •  L53/JG  -  x5  -  x18  (74 

T78  •  L37/JG  =  x18  -  x20  (75 

T45  *  L54/JG  "  x5  "  x19  (76 

T85  *  WJG  "  x19  "  x21  (77 

d.   Equations  Involving  Hooke's  Law 

The  tractor  tires  are  modeled  as  linear  springs  with  viscous  damping  as 
shown  in  Figure  8.  By  Hooke's  law,  the  force  at  the  tire  to  ground  contact 
point  of  the  i6"  mass  in  the  j  direction  is  given  by 
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Fij  =  Hj  '  (xij  "  yij)  +  Cij  '  (iij  "  yij)   i-1.2,3,4,7,8;3-l,2,3   (78 

Applying  (78)  to  the  six  tire  to  ground  contact  points  in  each  coordinate 
direction  yields 

Fll  =  ~(xl  +  (D53  +  D62)  •  x5  +  (Y  +  D99  +  D13)  •  x6  +  Du  •  x8  - 
yil)  •  KU  -(%  +  (D53  +  D62)  •  x5  +  (Y  +  D99  +  D13)  • 
x6  +  D14  •  x8  -  yu)  •  Cu  (79 

F12  =  -(x2  -  D53  •  x4  +  D51  •  x6  -  (D62  +  D14)  ■  x7  -  y12)  •  K12 

-(x2  -  D53  •  i4  +  D51  •  i6  -  (D62  +  D14)  •  i7  -  y12)  •  C12   (80 

F13  =  _(x3  "  Y  •  x4  -  D51  •  x5  -  (D99  +  D13)  •  x?  -  y13)  •  K13 

-(x3  -  Y  •  i4  -  D51  •  i5  -  (D99  +  D13)  •  x7  -  y13)  •  C13    (81 

F21  =  "(xl  +  (°53  +  D62>  '  x5  "  (°99  +  D22  "  *>  '  x6  +  D24  '  x9  " 
y2X)  •  K21  -(Xl  +  (D53  +  D62)  •  x5  -  (D99  +  D22  -  Y)  • 

x6  +  D24  -  x9  _  hi)    '  C21  (82 

F22  "  "(x2  "  D53  *  x4  +  D51  -  x6  "  (D62  +  D24>  '  x7  "  3^j)  '  K22 

-(x2  -  D53  ■  x4  +  D51  •  x6  -  (D62  +  D24)  •  x7  -  y22)  •  C22  (83 

F23  =  _(x3  "  Y  •  x4  -  D51  •  x5  +  (D99  +  D22)  ■  x7  -  y23)  •  K23 

-(x3  -  Y  •  x4  -  D51  •  x5  +  (D99  +  D22)  •  i7  -  y23)  •  C23    (84 

F31  =  "(x14  +  D34  '  x18  "  ?31>  '  K31  "  (x14  +  D34  *  x18  "  ?31>  '  C31   (85 
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F32  "  ~(x2  "  D57  '  x4  "  X  '  x6  "  D34  '  x26  "  y325  '  K32 

"(x2  "  D57  *  x4  "  X  •   x6  "  D34  '  x26  "  y32>  '  C32 


(86 


F33  -  -(x1Q  -  y33)  •  K33  -  (i1Q  -  y33)  •  C33 


(87 


F41  "  "(x15  +  D44  •  x19  "  y41>  *  K41  "  (*15  +  D44  '  x19  "  y41>  '   C41   (88 


F42  "  ~(x2  "  D57  '  x4  "  X  '  x6  "  D44  '  x27  "  y42>  '  K42 
-(x2  -  D57  •  x4  -  I  •  x6  -  D44  •  i27  -  y42)  • 


'42 


(89 


F43  =  -(xll  "  y43)  *  K43  "  (xll  "  y43>  * 


'43 


(90 


F71  =  "(x16  +  D74  *  x20  "  y71}  -  K71  "  (x16  +  D74  *  x20  "  y71>  "  C71   (91 


F72  =  "(x2  "  D57  •  x4  "  X  *  x6  "  D74  "  x28  "  y72>  '  K72 

-(x2  -  D57  •  i4  -  X  •  x6  -  D?4  •  x2g  -  y?2)  •  C?2 


(92 


F73  "  "(x12  "  y73>  •  K73  "  (x12  "  y73>  -  C73 


(93 


F81  =  "(x17  +  D84  "  x21  "  y81>  *  K81  "  (x17  +  D84  '  x21  "  y81>  '  C81   (94 


F82  '  "(x2  "  D57  *  x4  "  X  '  x6  "  D84  '  x29  "  y82)  '  K82 

-(x2  -  D57  •  x4  -  X  •  xg  -  Dg4  •  x29  -  yg2)  •  Cg2 


(95 


F83  =   ^x13  "  y83^  '  ^3  "  ^x13  "  y83^  ' 


'S3 


(96 
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e.  Reduction  of  the  Equations 

Equations  (4a)  through  (7f),  (9)  through  (72),  (74)  through  (77),  and 
(79)  through  (96)  constitute  one  hundred  ten  equations  in  one  hundred  ten 
unknowns.  Because  of  their  simplicity,  equations  (4a)  through  (7f)  are 
substituted  into  equations  (9)  through  (96)  where  applicable.  This  reduces 
the  number  of  equations  and  unknowns  to  eighty-six.  The  resulting  equations 
are  cast  into  matrix  form 


F.l 

*2 

*1 

*2 

h 

F 
86  x   57 

F81 

- 

A 
86  x   29 

k'3 

+ 

B 
86  x   29 

'] 

+ 

c 

86  x   29 

?3 

+ 

D 
86  x  39 

y18 

T33 

*2» 

i2, 

•« 

he 

FDBj 

FD»2 

FDB3 

(97 


Equation  (97)  is  rearranged  by  subtracting  A  from  either  side  and  com- 
bining B,  C  and  D  to  form 


FA 

F 

CBD 

X 

86  x 

36 

T 

X 

86  x    97 

X 

y 

y 

FDB 

(98 
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In  this  form,  the  eighty-six  unknowns  are  solved  for  using  the  Gauss- 
Jordan  reduction.  This  method  for  solving  simultaneous  equations  is  outlined 
in  reference  [3],  and  a  listing  of  the  computer  program  which  performs  the 
procedure  is  in  Appendix  A. 

After  reduction,  the  equations  are  of  the  form 


X 

F 

X 

T 

= 

CBD' 

y 

X 

y 

FDB 

(99 


where  I  is  the  identity  matrix  and  CBD1  is  a  recombination  of  CBD. 

The  last  twenty-nine  rows  of  (99)  are  the  differential  equations  of 
motion  for  the  system,  and  may  be  cast  into  the  form 


G 
19   x  57 


y 
y 

FDB 


(100 


where  G  is  the  last  twenty-nine  rows  of  CBD',  and  is  in  a  form  which  may  be 
integrated  numerically. 


B.  Derivation  Using  the  Energy  Method 

The  equations  of  motion  are  derived  by  the  energy  method,  to  be  used  as 
a  check  for  the  equations  found  in  part  A.  With  this  method,  expressions  for 
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the  kinetic,  potential,  and  dissipative  energy  are  formulated  and  used  to 
determine  the  equations  of  motion  with  the  Lagrange  equation  for  non- 
conservative  systems. 

As  in  part  A,  the  system  may  be  forced  by  either  a  drawbar  load  or  by 
the  surface  being  traversed.  The  drawbar  load  is  input  using  the  principle 
of  virtual  work,  and  the  surface  is  input  by  subtracting  the  profile  seen  by 
each  tire  in  each  coordinate  direction  from  the  applicable  displacement  in 
the  potential  energy  term. 

a.  The  Lagrange  Equation 

For  a  non-conservative  system  with  n  degrees  of  freedom,  the  Lagrange 
equation  takes  on  the  form 

d/dt  (3T/3xr)  -  3T/3xr  +  3F/3*r  +  3V/3xr  =  Qr     r-1,2,3  ...  n      (101 

where  T  is  kinetic  energy,  V  is  potential  energy,  and  F  is  the  dissipation 
function  due  to  damping  in  the  tires.   Qr  are  generalized  forces  acting  on 
the  system  that  are  not  included  in  V  or  F.  For  the  system  being  modeled 
here,  they  are  the  drawbar  forces. 

Application  of  equation  (101)  yields  a  set  of  n  second  order  differential 
equations.   Reference  [S]  outlines  the  formulation  of  the  Lagrange  equation 
and  its  application  to  both  conservative  and  non-conservative  systems. 

b.  Kinetic  Energy 

The  kinetic  energy  of  the  system  is  given  by 

T  =  h  •   I   Sj  •  x±2  (102) 
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t  -  H  '   &L   •    [(*!  +  (Y  +  D99  +  D13)  •  x6  +  (D53  +  Dfi2)  •  x5)2  +  C*j  + 
D51  *  *6  "  D53  •  H  ~   D62  '  *7>2  +  (*3  -  Y  •  x4  -  D31  •  i5  -  (D99  + 
D13)  •  x7)2]  +  Jlx  •  i72  +  Jn  •  x82  +  Jlz  •  x62  +  M2  •  [(Xl  +  (D53  + 
D62)  •  *5  -  (D99  +  D22  -  Y)  •  x6)2  +  C*2  +  D51  •  *6  "  D53  *  *4  "  D62  " 
x7)2  +  (x3  -  Y  •  x4  -  D51  •  *5  +  (D99  +  D22)  •  i7)2]  +  J2X  •  x?2  + 
J2y  •  i92  +  J2Z  •  i62  +  M3  •  [xu2  +  Cij  -  I  •  x6  -  D5?  •  x4)2  +  x102] 

+  J3X  '  *262  +  J3Y  '  *182  +  J3Z  '  *222  +  M4  '  t*152  +  (*2  "  *  •  *6  - 
D5?  •  x4)2  +  xu2]  +  J4X  •  x272  +  J4y  •  xig2  +  J4Z  •  x232  +  M5  •  [xx2 
+  x22  +  x32]  +  J5X  •  x42  +  J5Y  •  x52  +  J5Z  •  A62  +  Mg  •  [(X]_  +  (D53  + 

D63>  *  *5  +  *  *  *6>2  +  C*2  "  D53  •  *«  +  D51  *  *6  "  D63  '  *7>2  +  C*3  " 

Y  •  x4  -  D51  •  x5)2]  +  J6X  •  i72  +  J6Y  •  x52  +  Jgz  •  x62  +  M7  •  [x162 

+  (x2  -  D57  •  x4  -  X  •  x6)2  +  x122]  +  J7X  •  x282  +  J7Y  •  x202  +  J7Z  * 
*242  +  M8  *  ^172  +  C*2  "  D57  *  x4  -  X  •  x6)2  +  *132]  +  Jgx  •  *292  + 

J8Y  *  *212  +  J8Z  *  *252J  (103 

c.  Potential  Energy 

The  potential  energy  of  the  system  comes  from  two  sources,  the  deflection 
of  the  tires,  and  the  deflection  of  the  axles.  The  potential  energy  of  the 
deflected  tires  is  given  by 
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vx  =  h  •   I   %   xt2  (104 

Tt  -  H  '    [Ku  •  [*i  +  (D53  +  ^  ■  x5  +  (Y  +  D99  +  D13)  •  x6  +  D14  • 
x8  -  yu]2  +  Ku  •  [x2  -  D53  •  x4  +  D51  •  x6  -  (D62  +  Du)  •  x?  - 
y12]2  +  K13  •  [x3  -  Y  •  x4  -  D51  •  x5  -  (D99  +  Dn)  •  x7  -  y^]2 
+  K21  •  [xx  +  (D53  +  D62)  •  x5  -  (D99  +  D22  -  Y)  •  x6  +  D24  •  x9 
-  721]2  +  K22  •  [x2  -  D53  •  x4  +  D51  •  xfi  -  (Dg2  +  D24)  •  x?  -  y22]2 
+  K23  •  [x3  -  Y  •  x4  -  D51  •  x5  +  (D99  +  D22)  •  x7  -  y23]2  +  K31  • 
[x14  +  D34  •  X18  "  Y31^2  +  K32  *  £x2  "  D57  '  x4  "  *  '   x6  "  D34  ' 

x26  "  y32^  +  K33  •  ^x10  "  ?33]2  +  K41  '  Cx15  +  D44  -  x19  "  ^41^ 
+  K42  •  [x2  -  D57  •  x4  -  X  •  x6  -  D44  •  x27  -  y42]2  +  K43  •  [xu  - 

y43]2  +  K?1  •  [xlg  +  Dy4  •  x2Q  -  y?1]2  +  Ky2  •  [x2  -  D57  •  x4  -  X  • 
x6  "  D74  *  x28  "  772^  +  K73  '  tx12  "  y73^  +  K81  *  tx17  +  D84  ' 

x21  "  y81^  +  ^2  '  l-x2  "  D57  '  x4  "  X  '  x6  "  D84  '  x29  "  y82^ 

+  Kg3  •  [x13  -  y83]2]  (105 

A  discussion  of  the  potential  energy  in  the  deflected  axle  is  given  in 
Appendix  B.  This  energy  is  given  by 
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V2  =  (6EI/L533)  •  (-  x3  +  (D56  +  L53)  •  x4  -  X  •  x5  +  x10  +  D33  .  x26)2 
+  (6EI/L532)  •  (-  x3  +  (D56  +  L53)  •  x4  -  X  •  x5  +  x10  +  D33  •  x20)  • 
(-  x4  +  x26)  +  (2EI/L53)  .  (-  x4  +  x26)2  +  (6EI/L373)  •  (-  x10  + 
x12  +  (D32  +  L37)  •  x26  +  D73  •  x28)2  +  (6EI/L372)  •  (-  x10  +  x12 
+  (D32  +  L37)  •  x26  +  D73  .  x2g)  •  (-  x26  +  x28)  +  (2EI/L37)  •  (-  x26 
+  x2g)2  +  (6EI/L543)  •  (-  x3  -  (D55  +  L54)  •  x4  -  X  .  x5  +  xu  - 
D42  •  x27)2  -  (6EI/L542)  •  (-  x3  -  (D55  +  L54)  .  x4  -  X  .  x5  +  xu 
"  D42  •  x27>  '  C-  X4  +  x27)  +  (2EI/L54)  •  (-  x4  +  x27)2  +  (6EI/L843) 
•  C-  *U  +  X13  "  (^43  +  L84)  -  x27  -  D82  *  x29^2  "  (6EI/L842)  . 
(-  xu  +  x13  -  (D43  +  L84)  •  x27  -  D82  •  x29)  •  (-  x27  +  x29)  + 
(2EI/L84)  •  (-  x27  +  x29)2  +  (6EI/L533)  •  (-  xx  -  D57  •  x5  -  (D56  + 
L53)  .  x6  +  x14  -  D33  •  x22)2  -  (6EI/L532)  •  (-  *i   -  D57  •  x5  - 
(D56  +  L53)  •  x6  +  x14  -  D33  •  x22)  •  (-  x6  +  x22)  +  (2EI/L53)  • 
(-  x6  +  x22)2  +  (6EI/L373)  •  (-  x14  +  x16  -  (D32  +  L37)  •  x22  - 
D73  •  x24)2  "  (6EI/L372)  .  (-  x14  +  x16  -  (D32  +  L37)  .  x22  -  D73  • 
x24)  •  (-  x22  +  x24)  +  (2EI/L37)  .  (-  x22  +  x24)2  +  (6EI/L543)  • 
(-  x1   -   D57  .  x5  +  (D55  +  L54)  •  x6  +  x15  +  D42  •  x23)2  +  (6EI/L542)  • 
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(-  xx  -  D57  •  x5  +  (D55  +  L54)  •  x6  +  x15  +  D42  •  x23)  •  (-  x6  +  x23) 
+  (2EI/L54)  •  (-  x6  +  x23)2  +  (6EI/LgA3)  •  (-  x15  +  x17  +  (D43  +  Lg4)  • 
x23  +  D82  •  x25)2  +  (6EI/Lg42)  •  (-  x15  +  x17  +  (D43  +  Lg4)  •  x23  + 
D82  '  «23)  '  ("  x23  +  x25>  +  (2EI/L84)  •  (-  x23  +  x25)2  +  h  '    (JG/L53) 
•  C-  Xj  +  x18)2  +  h  •    (JG/L3?)  •  (-  xlg  +  x2Q)2  +  h   •  (JG/L59)  •  (-  x5 
+  x19)2  +  H  •    (JG/L84)  •  (-  x19  +  x21)2  (106 

The  total  potential  energy  Is  then 

VT  =  Vx  +  V2  (107 

d.  Dissipative  Energy 

The  dissipative  energy  due  to  the  damping  in  the  tires  is  given  by 

f  =  h  •   z:  c±  •  x±  (108 

F  -  h  '    [Cu  •  [xx  +  (D53  +  D62)  •  x5  +  (Y  +  D99  +  D13)  •  x6  +  D14  • 

*8  "  hl^2  +  C12  *  t*2  "  D53  '  *4  +  D51  -  46  "  <D62  +  D14>  *  *7  " 
y12]2  +  C13  •  [x3  -  Y  •  x4  -  D51  •  x5  -  (D99  +  D13)  •  x7  -  y13]2 
+  C21  •  [ix  +  (D53  +  D62)  •  i5  -  (D9g  +  D22  -  Y)  •  x6  +  D24  •  x,  - 
y21]2  +  C22  •  [x2  -  D53  •  x4  +  D51  •  x6  -  (Dg2  +  D24)  •  x?  -  y22]2 


+  C23  •  [x3  -  D51  •  x5  -  Y  •  x4  +  (D99  +  D22)  •  x7  -  y23]2  +  C 
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C*u  +  °34  -  *16  "  ^31]2  +  C32  '  [*2  "  D57  '  *4  "  *  '   *6  "  D34  '  *26 
-  *32]2  +  C33  •  [410  -  y33]2  +  c41  *  t*U  +   D44  *  *19  "  hl¥  * 

C42  -  t*2  "  D57  -  *4  "  *  "  *6  "  D44  '  *27  "  ^42^  +  C43  '  C*U  _ 
J43]2  +  C71  •  [*16  +  D?4  •  x2Q  -  yyl]2  +  C72  •  [*2  -  D57  •  *4  - 
X  •  46  -  D„  •  £28  -  y72]2  +  C73  •  [i12  -  y73]2  +  C81  •  [i17  + 

D84  •  *21  "  ^81^  +  C82  *  C*2  "  D57  '  *4  _  J  "  *6  "  D84  '  *29  " 

^8232  +  C83   •    [*13  -  f83J]  (109 

e.  Generalized  Forces 

The  generalized  forces  acting  on  the  system  are  due  to  the  components  of 
the  drawbar  load,  namely  FDB^,  FDB2  and  FDB3.  Using  the  principle  of  virtual 
work  as  outlined  in  reference  [8] 

6W  =  E  Q±  •  JJEj  (110 

=  FDB1-6x1  +  FDB,-(5x2  +  FDB3'<5x3  +  (-FDB2"DB3  -  FDB3"Y)'5x4 
+  (FDB^DBj  +  FDB3'DB1)-«x5  +  (FDB^T  -  FDB^DB^  -Sxg  (111 

Therefore 

Ql   =  FDB^  Q5  =  FDB1-DB3  +  FDB3-DB]_ 

Q2  =  FDB2  Qg  =■  FDB^Y  -  FDB2'DBX 

(112  a-g 
Q3  =  FDB3  Q7  through  Q29  =  0 

Q4  -  FDB2-DB3  -  FDB3-Y 
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f.  Reduction  of  the  Equations 

Substitution  of  the  expressions  for  T,  V,  F  and  QR  into  equation  (101) 
with  n  -  29  yield  the  equations  of  motion  for  the  system 

(M1  +  M2  +  M5  +  M6)  •  *j  +  [D53  •  (Mj_  +  M2  +  M6)   +  D62  •  (Mj_  +  M2) 

+  D63  •  M6]  •  i5  +  [?  •  (M]_  -  M2  +  Mg)  +  D99  •  (M]_  -  M2)  +  D13  •  Mj_ 

-  D22  •  M2]  •  *g  +  (CU  +  C21)  •  %  +  (D53  +  D62)  •  (Cu  +  C21)  •  x5 

+  [Cn  •  (Y  +  D99  +  D13)  +  C21  •  (Y  -  D99  -  D22)]  •  *g  +  Cjj,  ■  D14  ■  xg 
+  C21  •  D24  •  *g  -  Cn  •  yn  -  C21  •  y21  +  (Kn  +  K21  +  KAj_  +  KA3)  •  xL 
+  [(D53  +  D62)  •  (Kn  +  K21)  +  D57  •  (KAX  +  KA3)]  •  x5  +  [Ku  •  (Y  + 
D99  +  D13)  +  K21  •  (Y  -  D99  -  D22)  +  KAX  •  (D56  +  L53)  -  KA3  •  (D55  + 
L54)  -  KA5  +  KA7]  •  x6  +  Kn  •  D14  •  x8  +  K21  •  D24  •  xg  -  KAj_  •  x14 

-  KA3  •  x15  +  (KAX  •  D53  +  KA5)  ■  x22  +  (-KA3  •  D42  -  KA7)  •  x23  - 

Kll  *  ni  "  K21  *  ?21  =  roBl  (113 

(Mx  +  M2  +  M3  +  H4  +  M5  +  M6  +  My  +  M8)  •  x2  +  [-  D53  •  (M]_  +  M2  +  Mg) 

-  D57  •  (M3  +  M4  +  M7  +  M8)]  •  x4  +  [D51  •  (Mx  +  M2  +  Mg)  -  X  •  (M3  + 
M4  +  M7  +  Mg)]  •  X6  +  [-  D62  •  0*1  +  «2>  "  D63  "  •%]  '  X7  +  (C12  + 
C22  +  C32  +  C42  +  Cy2  +  C82)  •  x2  +  [-  D53  •  (C12  +  C22)  -  D57  •  (C32 
+  C42  +  C?2  +  C82)]  •  x4  +  [D51  •  (C12  +  C22)  -  X  •  (C32  +  C42  +  C,2 
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+  C82)]  •  x6  +  [-  C12  •  (D62  +  D14)  -  C22  •  (D62  +  D24)]  •  i?  -  C32  • 
D34  '  x26  "  c42  '  D44  '  x27  "  c72  '  D74  '  x28  "  c82  "  D84  '  x29  " 
C12  '  *12  "  C22  '  ?22  "  C32  '  ?32  "  C42  "  ?42  "  C72  '  ^72  "  C82  *  ^82 
+  (K^2  "*"  ^22  +  K32  +  K42  +  ^72  "**  ^82^  "  x2  "**  '-"  ^53  "  ^12  +  ^22'  ~~ 
D57  •  (K32  +  K42  +  K72  +  K82)]  •  x4  +  [D51  •  (K12  +  K22)  -  X  •  (K32 
+  K42  +  K72  +  K82)]  •  x6  +  [-  K12  •  (D62  +  D14)  -  K22  •  (D62  +  D24)] 

"  x7  "  K32  '  D34  "  x26  "  K42  '  D44  '  x27  "  K72  '  D74  '  x28  "  ^2  ' 

D84  *  x29  "  K12  *  ^12  -  K22  *  ?22  "  K32  "  ?32  "  K42  '  ?42  "  K72  *  ?72 

"  ^2  *  ?82  =  FDB2  (114 


(Mx  +  M2  +  M5  +  Mg)  •  x3  -  7  •  (Mj^  +  M2  +  >%)  •  »4  -  D51  •  (Mx  +  M2 
+  M6)  •  X5  +  [D99  •  (M2  -  Mx)  -  D13  •  MX  +  D22  •  M2]  •  »7  +  (C13  + 
C23)  •  i3  -  Y  •  (C13  +  C23)  •  x4  -  D51  •  (C13  +  C23)  •  x5  +  [-  C13  - 
(D99  +  D13)  +  C23  •  (D99  +  D22)]  •  x7  -  C13  •  y13  -  C23  •  y23  + 
(K13  +  K23  +  KAX  +  KA3)  •  x3  +  [-  Y  •  (K13  +  K23)  -  (Djg  +  L53)  • 
KAX  +  (D55  +  L54)  •  KA3  +  KA5  -  KAy]  •  x4  +  [-  D51  •  (K13  +  K23)  + 
X  •  (KAX  +  KA3)]  •  x5  +  [-  ^3  •  (D99  +  D13)  +  K^  •  (D^  +  0^)]  •  x? 
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-  KAX  •  x10  -  KA3  •  jcn  +  (-  KAX  •  D33  -  KA5)  •  x2g  +  (KA3  •  D42 

+  KA7)  •  x27  -  K13  •  y13  -  K23  •  y23  =  m3  (115 

[-  D53  •  CMj,  +  M2  +  Mg)  -  D5?  •  (M3  +  M4  +  M7  +  Mg)]  •  x2  -  T  •  (Mx 
+  M2  +  M6)  •  x3  +  [(D532  +  Y2)  •  (Mx  +  M2  +  Mg)  +  D572  •  (M3  +  M4  + 
M7  +  Mg)  +  J5X]  •  X4  +  Y  •  D51  •  (Mj_  +  M2  +  Mg)  ■  Xj  +  [-  Djj  •  D53  • 
(M]_  +  M2  +  Mg)  +  X  •  D57  •  (M3  +  M4  +  M?  +  Mg)]  •  Xg  +  [D53  •  Dg2  • 
(Mj_  +  M2)  +  D53  •  Dg3  •  Mg  +  Y  •  D99  •  (Mx  -  M2)  +  Y  •  Djj  •  Mj_  - 
Y  •  D22  •  M2]  •  x?  +  [-  D53  •  (C12  +  C22)  -  D5?  •  (C32  +  C42  +  C?2 
+  Cg2)]  •  x2  -  Y  •  (C13  +  C23)  •  x3  +  [D532  •  (C12  +  C22)  +  Y2  • 
(C13  +  C23)  +  D5?2  •  (C32  +  C42  +  C72  +  C82)]  •  x4  +  D51-   Y  •  (C13 
+  C23)  •  x5  +  [-  D51  •  D53  •  (C12  +  C22)  +  X  •  D57  •  (C32  +  C42  +  C72 
C82)]  •  *g  +  [D53  •  CCD62  +  Du)  •  C12  +  (Dg2  +  D24)  •  C22)  +  Y  • 
((D99  +  D13)  •  C13  -  (D99  +  D22)  •  C23)]  •  x7  +  C32  •  D57  •  D34  • 
x2g  +  C42  •  D57  •  D44  •  x27  +  C72  •  D57  •  D74  •  x2g  +  Cg2  •  D57  • 
D84  •  x29  +  C12  •  D53  •  y12  +  C13  •  Y  •  y13  +  C22  •  D53  •  y22  + 
C23  •  Y  •  y23  +  c32  •  D5?  ■  y32  +  C^    •   ^    •  ^   +  C?2  •  D„  •  yJ2 
+   C82  •  D57  •  y82  +  [-  D53  •  (K12  +  K22)  -  D57  •  (K32  +  K42  +  K?2  + 
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K82)]  •  x2  +  [-  Y  •  (K13  +  K23)  -  (D56  +  L53)  •  K*j  +  (D55  +  L54)  . 
KA3  +  KA5  -  KA?]  •  x3  +  [D532  •  (K12  +  K22)  +  Y2  ■  (Kjj  +  K23)  + 
D572  •  (K32  +  K42  +  Ky2  +  Kg2)  +  (D56  +  L53)2  •  KAj_  +  (Dgg  +  Lg4)2  • 
KA3  -  2  •  KA5  •  (D56  +  L53)  -  2  •  KAy  •  (Dgg  +  Lg4)  +  KAg  +  KAj^]  • 
x4  +  [D51  •  ?  •  (K13  +  K23)  -  X  •  ((D56  +  L53)  •  KAX  -  (D55  +  Lg4)  • 
KA3  -  KA5  +  KA7)]  •  x5  +  [-  D51  •  D53  •  (K12  +  K22)  +  X  •  Dg7  • 
(K32  +  K42  +  K72  +  Kg2)]  •  x6  +  [D53  •  ((D62  +  Du)  •  K12  +  (Dg2  + 
D24)  •  K22)  +  Y  •  ((Dg9  +  D13)  •  K13  -  (D99  +  D22)  •  K23)]  •  x?  + 
[KA-l  •  (D56  +  L53)  -  KA5]  •  x1Q  +  [-  KA3  •  (D55  +  L54)  +  KA?]  •  xu 
+  [D34  •  D57  •  K32  +  D33  •  (D56  +  L53)  •  KAj_  -  KAg  •  (D33  -  Dg6  - 
Lg3)  -  KA9]  •  x26  +  [D44  •  Dg7  •  K42  +  D42  •  (Dgg  +  Lg4)  •  KA3  + 
KA?  •  (-  D42  +  Dgg  +  Lg4)  -  KAjj]  •  x2 ?  +  D?4  •  Dg7  •  K?2  •  x28  + 

D84  '  D57  '  K82  '  x29  +  K12  '  D53  '  y12  +  K13  '  Y  '  y13  +  K22  '  D53  ' 
y22  +  K23  •  Y  •  y23  +  K32  •  Dg?  •  y32  +  K42  ■  Dgy  ■  y42  +  K?2  •  Dg?  • 

y72  +  K82  •  Dg7  •  y82  =  -DB3  •  FDB2  -  Y  •  FDB3 
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[(D53  +  D62)    •    (Mx  +  M2)   +   (D53  +  D63)    •   M6]    •   x1  -  D51    •    (f^  +  M2 
+  M6)    •    x3  +  D51    •  f   •    (Mx  +  M2  +  M6)    •   x4  +   [(D53  +  Dg2)2    •    (Mx  + 
M2)   +   (D53  +  D63)2    •   Mg  +  D512    •    (Mx  +  M2  +  Mg)   +  JJY  +  JfiY]    •   x5 
+  [(D53  +  D62)    •    (ff  +  D99  +  D13)    •   Hx  +   (Y  -  D99  -  D22)    •   M2)   + 
(D53  +  D63)    •   Y    •   Mg]    •   x6  +  D51    •    ((D99  +  D13)    •   Mx  -    (D99  +  D22)    • 
M2)    •   x7  +   (D53  +  D62)    •    (Cu  +  C21)    •   ix  -  D51    •    (C13  +  C23)    •   *3 
+  Y   •    D51    •    (C13  +  C23)    •   x4  +   [(D53  +  D52)2    •    (Cu  +  C21)   +  D512    • 
(C13  +  C23)]    •   *5  +   (C53  +  D62)    •    [(Y  +  D99  +  D13)    •    Cu  +   (Y  -  Dg9 
-  D22)    •    C21]    •   x6  +  D51    •    ((D99  +  D13)    •   C13   -    (D99  +  D22)    •    C23)    • 
x7  +  Cn    •    D14    •    (D53  +  D62)    •    ig  +  C21    •    D24    •    (D53  +  D62)    •   ig  - 

Cll    *    <D53  +  D62>    '   Hi  +  C13    •    D51    *    *13   "  C21    '    <D53  +  D62>    '    hi 
+  C23    •    D51    •    y23  +   [(D53  +  D62)    •    (Ku  +  K21)    +  D57    •    (KA]_  +  KA3)]    • 

xl  +   C-  D51    •    (K13  +  K23>   +  *   •    (KAl  +  KA3)]    *   x3  +  [d51    •   Y    •    (K13  + 
K23)    -  X    •    ((D56  +  L53)    •    KA-l   -    (D55  +  L54)    •   KA3  -  KA5  +  KA?)]    •   x4 
+  [(D53  +  D62)2    •    (Ku  +  K21)   +  D512    •    (K13  +  K23)   +   (X2  +  D572)    • 
(KAX  +  KA3)   +  KA13  +  KA15]    •   X5  +   [(D53  +  D62)    •    ((Y  +  Dgg  +  D13)    • 
Ku  +   (Y  -  D99  -  D22)    •    K21)   +  D5?    •    ((D56  +  L53)    •    KAT   -    (D55  +  L54) 
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KA3  -  KA5  +  KA?)]  •  x6  +  D51  •  «T>gg   +  D^)  •  K13  -  (Dgg  +  D22)  • 
K23)  •  x?  +  Du  •  (D53  +  D62)  •  Ku  •  xg  +  D24  •  (D53  +  Dfi2)  .  K^ 


Xq  —  X  •  KA*  ■  xiq  "  X  •  KAq  *  xll  ""  57  *  ^"i  *  x14  ~  57  *  *^"3 


x15  "  ^13  *  x18  "  ^15  "  x19  +  D57  "  (KA1  '  D33  +  ^S5  '  x22  + 


D5?  •  (-  KA?  -  KA3  •  D42)  •  x23  +  (-  X  •  D33  •  KAX  -  X  •  KAj)  •  x2fi 


(X  •  D42  •  KA3  +  X  •  KA7)  •  x27  -  K1]_  •  (D53  +  Dg2>  •  yn 


"13 


D51  •  yU  "  K21  •  <D53  +  D62>  •  *21  +  K23  *  D51  '  ?23  =  DB3  '  ™1 


+  DBj^  •  FDB3 


(117 


[(?  +  D99  +  D13)  •  Ml   +  (Y  -  D99  -  D22)  •  M2  +  Y  •  Mg]  •  Xl  + 
[D51  •  0*L  +  M2  +  Mg)  -  X  •  (M3  +  M4  +  M7  +  M8)]  •  x2  +  [-  D51  • 
D53  •  (J^  +  M2)  -  D51  •  D53  •  Mg  +  D57  •  I  •  (M3  +  M4  +  M7  +  Mg)]  • 
x4  +  [CD53  +  Dg2>  •  ((Y  +  D99  +  D13)  •  MX  +  (Y  -  D99  -  D22)  ■  M2) 
+  Y  •  (D53  +  Dg3)  •  Mg]  •  S5  +  [(Y  +  D99  +  D13)2  •  Mj.  +  (Y  -  D99  - 
D22)2  •  M2  +  D512  •  (Mx  +  M2  +  Mg)  +  X2  •  (M3  +  M4  +  M7  +  Mg)  +  Y2  • 

M6  +  J1Z  +  J2Z  +  J5Z  +  J6Z^  *  x6  +  ["  D51  '  D62  '  <M1  +  M2>  "  D51  ' 
Dg3  •  Mg]  •  X7  +  [CU  •  (Y  +  D99  +  D13)  +  C21  •  (Y  -  D99  -  D22)]  •  *j 

+  [D51  •  (C12  +  C22)  -  X  •  (C32  +  C42  +  C72  +  C82)]  •  *2  +  [-  D51  • 
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D53  *  (c12  +  C22)  +  X  •  D57  •  (C32  +  C42  +  C72  +  C82)]  ■  x4  +  (D53  + 
D62)  •  (Cu  •  (f  +  D99  +  D13)  +  C21  •  (Y  -  D99  -  D22))  •  x5  +  [Cu  • 
(Y  +  D99  +  D13)2  +  C21  •  (Y  -  D99  -  D22)2  +  D512  •  (C12  +  C22)  + 
X2  •  (C32  +  C42  +  C72  +  C82)]  •  x6  -  D51  •  ((D62  +  D14)  •  C12  + 
(D62  +  D24)  •  C22)  ■  i7  +  C1X  •  D14  •  (Y  +  D99  +  D13)  •  xg  +  C21  • 
D24  •  (Y  -  D99  -  D22)  •  x9  +  C32  •  X  •  D34  •  x26  +  C42  •  X  •  D44  • 
i27  +  C72  •  X  •  D74  •  x28  +  C82  •  X  •  D84  •  x29  -  Cu  •    (Y  +  D99  + 
D13)  •  yu  -  C12  •  D51  •  y12  -  C21  •  (Y  -  D99  -  D22)  •  y21  -  C22  • 
D51  *  ?22  +  C32  •  X  •  y32  +  C42  •  X  •  y42  +  C72  ■  X  •  y72  +  C82  • 
X  •  y82  +  [Ku  •  (Y  +  D99  +  D13)  +  K21  •  (Y  -  D99  -  D22)  +  KAX  • 
(D56  +  L53)  -  KA3  •  (D55  +  L54)  -  KA5  +  KA7]  •  xx  +  [D51  •  (K12  + 
K22)  -  X  •  (K32  +  842  +  K72  +  K82)]  •  x2  +  [-  D51  •  D53  •  (K12  +  K22) 
+  D57  •  X  •  (K32  +  K42  +  K72  +  K82)]  •  x4  +  [(D53  +  D62)  •  ((Y  +  Dgg 
+  D13)  •  Ku  +  (Y  -  Dg9  -  D22)  •  K21)  +  D57  •  ((D56  +  L53)  •  KA.J  - 
(D55  +  L54)  •  KA3  -  KA5  +  KA7)]  •  x5  +  [  ( Y  +  D99  +  D13)2  •  Ku  + 
(Y  -  D99  -  D22)2  •  K21  +  D512  •  (K12  +  K22)  +  X2  •  (K32  +  K42  +  K72 
K82)  +  ^D56  +  L53>2  -  K*!  +  (D55  +  L54>2  •  KA3  -  2  •  KA5  •  (D56  + 
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L53)  -  2  •  KA7  •  (D55  +  L54)  +  KA9  +  KAn]  •  x6  +  [-  D51  •  (D62  + 

D14>  -  K12  "  D51  "  (D62  +  D24>  '  K223  "  x7  +  D14  '  ff  +  D99  +  Dl3)  " 
Kll  '  x8  +  D24  •  (Y  -  D99  -  D22)  •  K21  •  x9  +  [-  KAX  •  (D56  +  L53) 
+  KA5]  •  x14  +  [KA3  •  (D55  +  L54)  -  KA7]  •  x15  +  [D33  •  (D56  +  L53) 

•  KAX  +  KA5  •  (-  D33  +  D55  +  L53)  -  KA9]  •  x22  +  [D42  •  (D55  +  L54) 

•  KA3  +  KA7  •  (-  D42  +  D56  +  L54)  -  KAn]  •  x23  +  X  •  D34  •  K32  •  x26 
+  X  •  D44  •  K42  •  x27  +  X  •  D74  •  K72  •  x28  +  X  •  D84  •  K82  •  x29  - 
Ku  •  (Y  +  D99  +  D13)  •  yu  -  K12  •  D51  •  y12  -  K21  •  (?  -  D99  -  D22) 

•  Y21  "  K22  "  D51  '  Y22  +  K32  •  X  •  y32  +  K42  •  X  •  y42  +  K72  •  X  • 

y72  +  K82  •  X  •  y82  =  7  •  YDB1   -  DBj^  •  FDB2  (118 


[-  D62  •  (M]i  +  M2)  -  D63  •  Mg]  •  x2  +  (D99  +  D13)  •  (M2  -  Mx)  •  X3 
+  tD53  *  D62  *  (Ml  +  M2)  +  D63  *  D53  •  M6  +  ^  •  (D99  +  D13)  •  (Mj.  - 
M2)]  •  x4  +  (D99  +  D13)  •  D51  •  (MX  -  M2)  •  x5  +  [-  D62  •  D51  •  (ML 
+  M2)  -  D63  •  D51  •  Kg]  •  x6  +  [D622  .  (MX  +  M2)  +  D632  .  M6  +  (D99 
+  D22)2  •  (Mi  +  M2)  +  Jlx  +  J2X  +  J6X]  .  x7  +  [-  C12  •  (D62  +  D14)  - 
C22  •  (D62  +  D24)]  •  x2  +  [-  C13  •  (D99  +  D13)  +  C23  •  (D99  +  D22)] 
i3  +  [D53  •  (C12  •  (D62  +  D14)  +  C22  •  (D62  +  D24))  +  T  •  (C13  (D99 


45 

+  D13)  -  C23  •  (D99  +  D22))]  •  x4  +  D51  •  (C13  •  (D99  +  D13)  -  C23  • 
(D99  +  D22))  •  i5   -  D51  •  (C12  •  (D62  +  D14)  +  C22  •  (Dg2  +  D24))  • 

h   +  Cc12  •  <D62  +  D14)2  +  c13  *  <D99  +  D13)2  +  c22  *  <D62  +  D24>2  ♦ 
C23  •  (D99  +  D22)2]  •  x7  +  C12  •  (D62  +  D14)  •  y12  +  C13  •  (D99  +  D13) 

"  ?13  +  c22  *  (°62  +  D24>  '  ?22  "  c23  '  <D99  +  D22>  '  ?23  +  i"   <D62  + 
D14)  •  K12  -  (D62  +  D24)  •  K22]  •  x2  +  [-  (D9g  +  D13)  •  K13  +  (D99  + 
D22)  •  K23]  •  x3  +  [D53  •  ((D62  +  D14)  •  K12  +  (D62  +  D24)  •  K22)  + 
T  •  ((D99  +  D13)  •  K13  -  (D99  +  D22)  •  K23)]  •  x4  +  D51  •  ((D99  + 
D13)  •  K13  -  (Dg9  +  D22)  •  K23)  •  x5  -  D51  •  ((D62  +  D14)  •  K12  + 
(D62  +  D24)  •  K22)  •  x6  +  [(D62  +  D14)2  •  K12  +  (D99  +  D13)2  •  K13  + 
(D62  +  D24)2  •  K22  +  (D99  +  D22)2  •  K23]  •  x7  +  K12  •  (D62  +  D14)  • 
y12  +  K13  •  (D99  +  D13)  •  y13  +  K22  •  (D62  +  D24)  •  y22  -  K23  •  (D99 
+  D22)  •  y23  =  0  (119 


J1Y  •  x8  +  Cll  •  D14  '  *1  +  Cll  '  D14  '  <D53  +  D62>  '  *5  +  CU  ' 
D14  *  «   +  D99  +  D13>  •  *6  +  Cll  *  D142  '  *8  '  Cll  '  D14  '  Hi  + 


hi   ■   D14  •  xl  +  hi    •  D14  *  (D53  +  D62>  '  x5  +  Kll  '  D14  '  (Y  + 


D99  +  D13>  •  x6  +  Kll  •  D14 


Kn  '  Du  •  rn  -  o 


(120 
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J2Y  •  x,  +  C21  •  Dw  •  xx  +  C21  •  D24  •  (D53  +  D62)  •  x5  +  C21 


D24  *  (Y  "  D99  "  D22>  '  *6  +  C21  "  D24 


x„  -  C 


9    21   "24   '21 


*»i  + 


K21  *  D24  *  xl  +  K21  *  D24  *  (D53  +  D62>  "  x5  +  K21  '  D24  '  (Y 


D99  "  D22)  *  x6  +  K21  •  D24 


x9  "  K21  '  D24  "  y21 


=  0 


M3  •  X10  +  C33  •  i10  -  C33  •  ?33  -  KAL  •  x3  +  [^  •  (D56  +  L53) 
KA5]  •  X4  -  M.J  •  X  •  x5  +  (K33  +  KAX  +  KA2)  •  x10  -  KA2  •  x12  + 
[KAX  •  D33  -  KA2  •  (D32  +  L37)  +  KA5  +  KAg]  •  x26  +  (-  KA2  •  D?3 
KAg)  •  x28  -  K33  •  y33  =  0 


(121 


(122 


M4  •  xn  +  C43  •  xu  -  C43  •  fi3  -  KA3  •  x3  +  [-  KA3  •  (D55  +  L54) 
+  KA7]  •  x4  -  KA3  •  X  •  x5  +  (K43  +  KA3  +  KA4)  •  xu  -  KA4  •  x13  + 
[-  KA3  •  D42  +  KA4  •  (D43  +  L84)  -  KA?  -  KAg]  •  x27  +  (KA4  •  D82  + 
KA8)  •  x29  -  K^3    •   y43  =  0  (123 


M7  •  x12  +  C73  •  x12  -  C73  •  y?3  -  KA2  •  x1Q  +  (K?3  +  KA2)  •  x12  + 

[KA2  •  (D32  +  L37)  -  KA6]  •  x26  +  (KA2  •  D73  +  KA6)  •  x28  -  K73  • 

773  ■  0  (124 
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M8  *  x13  +  c83  '  *13  "  C83  "  783  ~  ^4  '  xll  +  (K83  +  ^M~>    '   x13  + 

[-  KA4  •  (D43  +  L84)  +  KAg]  •  x27  +  (-  KA4  •  D82  -  KAg)  •  x29  - 

Kg3  •  y83  =  0  (125 

M3  •  «14  +  C31  •  %  +  C31  •  D34  •  x18  -  C31  •  y31  -  KAX  •  x1  - 

KA],  •  D57  ■  xj  +  [-  Kft  •  (D56  +  L53)  +  KA5]  •  Xg  +  (K31  +  KAX  + 


KA2)  •  xu  -  KA2  •  xlfi  +  K31  •  D34 


xlg  +  [-  KAj_  •  D33  +  KA2 


(D32  +  L37)  -  KA5  -  KAg]  •  x22  +  (KA2  •  D?3  +  KAg)  •  x24  -  K31  • 

731  =  °  (126 

M4  •  X15  +  C41  •  x15  +  C41  •  D44  •  x19  -  C41  •  y41  -  KA3  •  Xj  - 

KA3  •  D57  •  x5  +  [KA3  •  (D55  +  L54)  -  KA7]  •  x6  +  (K^  +  KA3  + 

KA4)  •  x15  -  KA4  •  x17  +  K41  •  D44  •  x19  +  [KA3  •  D42  -  KA4  • 

(D43  +  L84)  +  KA7  +  KA8]  •  x23  +  (-  KA4  •  Dg2  -  KA8)  •  x25  -  K41  • 

741  ■  0  (127 

M7  *  x16  +  c71  *  x16  +  c71  *  D74  '  x20  "  C71  '  *71  "  ^2  '  x14  + 

(K71  +  KA2)  •  x16  +  K71  •  Dy4  •  x20  +  [-  KA2  •  (D32  +  L37)  +  KAg] 

•  x22  +  (-  KA2  •  D73  -  KAg)  •  x24  -  K71  •  y71  =  0  (128 
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M8  •  x17  +  C81  •  x17  +  C81  •  D84  •  *21  "  C8l  •  y81  -  K^  •  x15  + 
(K81  +  KA4)  •  x17  +  K81  •  D84  •  x21  +  [KA4  •  (D43  +  L84)  +  KA8]  • 
x23  +  (KA4  •  D82  +  KA8)  •  x25  -  K81  •  y81  =  0  (129 

J3Y  •  x18  +  c31  *  D34  -  *14  +  c31  '  D342  '  *18  "  c31  *  D34  '  *31  " 

KA13  •  x5  +  K31  •  D34  •  x14  +  (K31  •  D342  +  KA13  +  KA14)  •  x18  - 

KA14  •  x20  -  K31  •  D34  •  y31  =  0  (130 

J4Y  •  X19  +  C41  *  D44  '  *15  +  C41  '  D442  '  *19  "  C41  '  D44  '  hi   " 
KA15  •  x5  +  K41  •  D44  •  x15  +  (K41  •  DA42  +  KA^  +  KA^)  .  ^  - 

**16  *  X21  "  K41  '  °44  •  ^1  =  °  <131 


J7Y  •  x2Q  +  Cn  •  D74  •  x16  +  Cn  •  D742  •  x2Q  -  C?1  •  Dy4  •  y?1  + 
K71  •  D74  *  x16  "  ^14  •  x18  +  <K71  *  D742  +  "lA5  "  x20  "  K71  '  D74 


y71 


(132 


J8Y  •  x21  +  Cgl  •  Dg4  •  x17  +  C81  •  D84   •  x21  -  C81  •  Dg4  •  y81  + 

K81  '  D84  •  x17  ~  **X6  *  x19  +  (K81  -  D842  +  "lS5  '  x21  "  %  '  D84 

*  y81  "  0  (133 
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J3Z  *  x22  +  0^1  *  D33  +  ^5)  *  xl  +  CK*i  •  D33  •  D57  +  D57  •  KA5) 

•  x5  +  [(D56  +  L53)  (KAX  •  D33  +  KA5)  -  KA5  •  D33  -  KAg]  •  x6  + 

[-  KAX  •  D33  +  KA2  •  (D32  +  L37)  -  KA5  -  KA6]  •  x14  +  [-  KA2  • 

(D32  +  L37)  +  KAg]  •  x16  +  [KAj_  •  D332  +  KA2  •  (D32  +  L37)2  + 

2  •  KA5  •  D33  -  2  •  KAg  •  (D32  +  L37)  +  KA9  +  KA10]  •  x22  +  [KA2  • 

D73  •  (D32  +  L37)  +  KAg  •  (-  D73  +  D32  +  L37)  -  KA10]  •  x24  =  0      (134 


J4Z  *  x23  +  (-  KA3  •  D42  -  KA7)  •  Xl  +  (-  KA3  •  D42  •  D57  -  D$7  • 
KA7)  *  x5  +  [(D55  +  L54)  •  (KA3  ■  D42  +  KA7)  -  KA?  ■  D42  -  KAn]  • 
x6  +  [KA3  •  D42  -  KA4  •  (D43  +  L84)  +  KA7  +  KAg]  •  x15  +  [KA^  • 
(D43  +  L84)  -  KA8]  •  x17  +  [KA3  •  D422  +  KA4  •  (D43  +  L84)2  + 
2  •  KA7  •  D42  -  2  •  KAg  •  (D43  +  L84)  +  KAU  +  KA12]  ■  x23  + 
[KA4  •  D82  •  (D43  +  L84)  -  KA8  •  (D82  -  D43  -  L84)  -  KA12]  •  x25 


=  0 


(135 


J7Z  '  S24  +  (^2  *  D73  +  "6>  "  x14  +  (-  KA2  •  Dy3  -  KAfi)  •  x16  + 

[(D32  +  L37)  •  (KA2  •  D73  +  KAg)  -  KAg  •  D73  -  KA10]  •  x22  +  (KA2  • 

D732  +  2  •  KAg  ■  D73  +  KA10)  •  x24  =  0  (136 


J8Z  *  x25  +  C-  KA4  •  D82  -  KA8)  •  x15  +  (KA4  •  D82  +  KAg)  •  x17  + 

[(D43  +  L84)  •  (KA4  •  D82  +  KA8)  -  KA8  •  D82  -  KA12]  ■  x23  +  (KA4  • 

D822  +  2  •  KAg  •  D82  +  KA12)  •  x25  =  0  (137 

J3X  *  x26  "  c32  *  D34  '  *2  +  c32  '  D34  '  D57  '  *4  +   c32  '  D34  '  '  * 

x6  +  C32  •  D342  •  A26  +  C32  •  D34  •  y32  -  K32  •  D34  •  x2  +  (-  KAX  • 

D33  -  KA5)  •  x3  +  [K32  •  D57  •  D34  +  (D56  +  L53)  •  (KAj.  •  D33  +  KA5) 

-  KA5  ■  D33  -  KAg]  •  x4  +  (-  KAX  •  D33  •  I  -  X  •  KA5)  •  x5  +  K32  • 

D34  •  X  •  x6  +  [KAX  •  D33  -  KA2  •  (D32  +  L37)  +  KA5  +  KAg]  •  x10  + 

[KA2  •  (D32  +  L37)  -  KAg]  •  x12  +  [K32  •  D342  +  KA]_  •  D332  +  KA2  • 

(D32  +  L37)2  +  2  •  KA5  •  D33  -  2  •  KAg  •  (D32  +  L37)  +  KA9  +  KA10]  • 

x2g  +  [KA2  •  (D32  +  L37)  •  D73  -  KAg  •  (D73  -  D32  -  L37)  -  KA10]  • 

x28  +  K32  •  D34  •  y32  =  0  (138 

J4X   •    x27   -  C42    •    D44    •    i2  +  C42    •    D44    •    D57    •    x4  +  C42    •    D44    •    X   • 

*6  +  C42    •    D442    '   i27  +  C42    '    D44    #    *42   "  ^2    '    D44    '    x2  +   ^3    " 
D42  +  KA;)    •   x3  +  [K42    .    D44    •    DS7   +   (D55  +  L54)    •    (KA3    •    D42  +  KAy) 


-  KA7 


D42   "  %]    *    x4  +   (^3    •    D42    •    X  +  X   •    KA7)    •    x5  +  K, 


42    - 
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X  •  x6  +  [-  KA3  •  D42  +  KA4  •  (D43  +  Lg4)  -  KA-,  -  KAg]  ■  xu  + 


[-  KA4  •  (D43  +  L84)  +  KAg]  •  x13  +  [K42  •  D442  +  KA3  •  D422  +  KA^    • 
(D43  +  L84)2  +  2  •  KA7  •  D42  -  2  •  KA8  •  (D43  +  L84)  +  KAU  +  KA12]  • 
x27  +  [KA4  •  (D43  +  L84)  •  D82  -  KAg  •  (D82  -  D43  -  Lg4)  -  KA12]  • 
x29  +  K^  •  D44  •  y42  =  0  (139 

J7X  '  x28  "  c72  •  D74  '  *2  +  c72  #  D74  '  D57  '  *4  +  c72  '  D74  *  *  ' 
*6  +  c72  •  D742  •  *28  +  c72  '  D74  "  ?72  "  K72  '   D74  *  x2  +  K72  '  D74  ' 
D57  "  x4  +  K72  •  D74  •  X  •  x6  +  (-  KA2  •  D73  -  KAg)  •  x10  +  (KA2  • 
D73  +  KAg)  •  x12  +  [KA2  •  D73  •  (D32  +  L37)  +  KAg  •  (D32  +  L37  -  D73) 

-  KA10]  '  x26  +  (K72  *  D742  +  ^2  '  D732  +  2  •  KAg  •  D73  +  KA10)  • 

x28  +  K72  •  D74  •  y72  =  0  (140 

J8X  -  x29  "  c82  *  D84  '  x2  +  c82  '  D84  '  D57  "  x4  +  c82  '  D84  -  X  • 
ig  +  C82  •  D842  •  i29  +  Cg2  •  Dg4  •  yg2  -  Kg2  •  D84  •  x2  +  K82  •  D84  • 
D57  •  x4  +  K82  •  Dg4  •  X  •  xg  +  (KA4  •  D82  +  KAg)  •  xu  +  (-  KA4  • 
Dg2  -  KAg)  •  x13  +  [KA4  •  Dg2  •  (D43  +  L84)  +  KA8  •  (D43  +  Lg4  -  Dg2) 

-  KA12]  •  x27  +  (Kg2  •  Dg42  +  KA4  •  D822  +  2  •  KAg  •  Dg2  +  KA12)  • 

x29  +  K82  -  D84  •  Y82  =  °  (141 
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where 


KAX  «  12EI/L533 
KA2  =  I2EI/L373 
KA3  -  I2EI/L543 
KA4  =  12EI/L843 
KA5  -  6EI/L532 
KAg  -  6EI/L372 
KA7  -  6EI/L542 
KA8  =  6EI/L842 


KA9  =  4EI/L53 
KA10  =  4EI/L37 
KAU  -  4EI/L54 
KA12  =  4EI/L84 
KA13  =  JG/L53 
KA14  =  JG/L37 
KA15  -  JG/L54 
KA16  =  JG/L84 


(142  a-p 


Equations  (113)  through  (141)  may  be  cast  into  the  matrix  foi 


X 

+ 

B 

X 

♦H 

X 

.[.] 

y 

y 
FDB 

(143 


Equation  (143)  may  be  rearranged  to  the  form 


CBD 


x 

y 
y 

FDB 


(144 


This  equation  is  now  in  a  form  which  allows  the  Gauss-Jordan  method  to 
be  used.   Upon  reduction,  equation  (144)  becomes 


CBD '  I 


y 
y 

FDB 


(145 


The  matrix  CBD'  is  identially  the  matrix  G  in  equation  (101) ,  allowing 
equation  (145)  to  be  rewritten  as 


X 

- 

G 

X 

X 

y 
y 

FDB 

(146 


VI.   SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  equations  of  motion  are  to  be  integrated  numerically  using  the  fourth 
order  Runge-Kutta-Gill  method.  This  method  is  used  to  integrate  first  order 
differential  equations  with  known  initial  conditions.  The  equations  of  motion 
for  this  system  are  second  order  differential  equations  which  may  be  trans- 
formed into  a  set  of  fifty-eight  first  order  equations  by  the  transformation 


*i  '  xi+29 


1  -  1,2,3 


29 
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Equations  (147)  and  (100)  or  (146)  constitute  the  set  of  fifty-eight 
equations  which  are  combined  in  the  single  matrix  equation 


58  x  58 


v;s 


If  l     W 

9J[29  x  29|  [: 


0    II    I 
29  x  29 1 | 29  x  29| 

G 
9  x  97 


1 

xi 

0 

29  x   29 

X58 

; 

yl 
FDB3 
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This  equation  is  now  in  a  form  which  may  be  integrated  to  give  the  dis- 
placement and  velocity  of  each  of  the  twenty-nine  x's  as  functions  of  time. 
Details  of  the  Runge-Kutta-Gill  numerical  integration  process  are  outlined 
in  reference  [10]. 
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It  should  be  noted  that  a  vector  x0  is  needed  to  start  the  Runge-Kutta- 
Gill  procedure.  The  first  twenty-nine  elements  of  this  vector  are  initial 
displacements,  and  the  last  twenty-nine  are  initial  velocities. 


VII.   THE  COMPUTER  PROGRAM 

The  computer  program  to  simulate  the  motion  of  a  tractor  is  broken 
into  two  parts.  The  first  part  compiles  the  various  tractor  parameters, 
formulates  the  equations  that  govern  the  system,  and  reduces  them  to  the 
twenty-nine  differential  equations  of  motion.  The  second  part  integrates 
the  equations  derived  in  Part  I, 

A.  Part  I 

Programs  to  derive  the  equations  of  motion  are  developed  using  both 
the  method  of  section  VA,  and  the  method  of  section  VB.  Both  programs  use 
the  same  tractor  parameters,  and  input  them  in  the  same  manner.  These 
parameters  are  used  to  evaluate  the  coefficients  multiplying  each  variable. 

A  list  of  the  parameters  needed,  and  the  order  in  which  they  are  read 
into  the  computer  is  shown  in  Table  I.   It  is  important  to  note  that  all 
distances  are  positive  according  to  Figure  6.  For  example,  the  center  of 
gravity  of  each  wheel  is  assumed  to  be  to  the  left  of  the  geometric  center 
of  the  wheel.   If  the  center  of  gravity  is  actually  to  the  right  of  the 
geometric  center,  that  distance  would  be  input  as  a  negative  value. 
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TABLE  I 


ORDER  IN  WHICH  THE  TRACTOR  PARAMETERS  ARE 
INPUT  TO  THE  COMPUTER 


Card Parameters 


1  "l*  J1X'  J1Y'  J1Z*  D13'  D14 

Kj,  J2X,  J2Y,  J2Z,  D22,  D24 

M3*  J3X'  J3Y'  J3Z'  D31*  D32*  D33'  D34 

4  M4,  J4X,  J4Y,  J4Z,  D41>  D42,  D43>  D44 

«%•  J5X,  J5Y.  J5Z.  D5i>  D53,  D55,  D56 
6                 D5? 

M6'  J6X'  J6Y'  J6Z'  D62>  D63'  D64'  D99 

M7,  J?x>  J7Y,  Jyz,  D71,  D73,  D74 

M8'  J8X'  J8Y'  J8Z'  D81>  D82*  D84 


9 


L37'  L53>  L54'  L84 

11  EI*,  JG* 

12  X,  Y,  Z,  DBj^,  DB3 

13  Kj^,  K12,  K13,  K^,  K22,  K23,  Zy±,   K32 

14  ^3'  K41'  K42'  K43'  K71'  K72'  K73'  % 

15  Kg2,  Kg3 

16  Cn>  C12,  C13,  C21,  C22,  C23>  C31,  C32 

17  C33*  C41«  C42'  C43'  C71"  C72*  C73«  C81 

18  C82,  Cg3 


*  These  values  are  input  in  E  format;  all  others  are  F  format. 
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a.  A  Program  Using  the  Method  of  Section  VA 

In  this  program,  the  coefficients  are  arranged  in  five  matrices  as  in 
equation  (97) .  These  five  matrices  are  cast  into  a  single  matrix  E  where 


HhlHHHl 


(149 


E  is  a  coefficient  matrix  which  multiplies  a  vector  consisting  of  the 
variables,  as  shown  in  Table  II. 

The  E  matrix  is  now  passed  into  the  Gauss-Jordan  subroutine  which 
reduces  it  to  the  equations  of  motion.   The  last  one  hundred  twenty-six 
columns  of  the  last  twenty-nine  rows  of  the  reduced  E  matrix  are  the  actual 
equations  of  motion.  These  rows  and  columns  are  put  into  a  matrix  G  which 
may  be  printed  out,  or  stored  on  magnetic  disk  for  use  in  the  integration 
program.  A  listing  of  this  program,  labeled  TRA,  is  in  Appendix  A. 

b.  A  Program  Using  the  Method  of  Section  VB 
This  program  arranges  the  coefficients  in  four  matrices  as  in 
equation  (143) .  These  matrices  are  cast  into  a  matrix  G  where 


[HHMH1 


(150 


This  matrix  is  passed  to  the  Gauss-Jordan  subroutine  for  reduction.  No 
further  manipulation  of  this  matrix  is  necessary,  as  it  contains  only  the 
equations  of  motion.  This  G  matrix  is  identical  to  that  of  TRA,  and  may 
also  be  printed  out  or  stored  on  disk  for  later  use.  A  listing  of  this 
program,  labeled  TRB,  is  also  in  Appendix  A. 
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TABLE  II 
ORDER  OF  THE  VARIABLES  IN  THE  VECTOR  MULTIPLYING  E 


Element 

Variable 

Element 

Variable 

Element 

Variable 

1 

Fll 

23 

F42 

45 

T35 

2 

»12 

24 

F43 

46 

T36 

3 

F13 

25 

F44 

47 

T37 

4 

F17 

26 

F45 

48 

T38 

5 

F18 

27 

F46 

49 

T39 

6 

F19 

28 

F47 

50 

T44 

7 

F21 

29 

F48 

51 

T45 

8 

F22 

30 

F49 

52 

T46 

9 

F23 

31 

F51 

53 

^47 

10 

F24 

32 

F52 

54 

T48 

11 

F25 

33 

F53 

55 

T49 

12 

F26 

34 

F71 

56 

T52 

13 

F31 

35 

F72 

57 

T53 

14 

F32 

36 

F73 

58 

*1 

15 

F33 

37 

F31 

59 

*2 

16 

F34 

33 

F82 

60 

*3 

17 

F35 

39 

F83 

61 

x4 

18 

F36 

40 

T17 

62 

*5 

19 

F37 

41 

H9 

63 

S6 

20 

F38 

42 

?24 

64 

X7 

21 

F39 

42 

T26 

65 

*6 

22 

F41 

44 

T34 

66 

x9 
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TABLE  II   (Continued) 


Element 

Variable 

Element 

Variable 

Element 

Variable 

67 

x10 

89 

x3 

111 

x25 

68 

*11 

90 

x4 

112 

x26 

69 

*12 

91 

x5 

113 

x27 

70 

x13 

92 

x6 

114 

x28 

71 

*14 

93 

x7 

115 

x29 

72 

x15 

94 

x8 

116 

xl 

73 

x16 

95 

x9 

117 

x2 

74 

'x17 

96 

x10 

118 

x3 

75 

x18 

97 

xll 

119 

x4 

76 

x19 

98 

x12 

120 

x5 

77 

x20 

99 

x13 

121 

x6 

78 

*u 

100 

x14 

122 

x7 

79 

K22 

101 

x15 

123 

x8 

80 

jc23 

102 

x16 

124 

x9 

81 

x24 

103 

x17 

125 

x10 

82 

x25 

104 

x18 

126 

xll 

83 

*26 

105 

x19 

127 

x12 

84 

*27 

106 

x20 

128 

x13 

85 

s28 

107 

x21 

129 

x14 

86 

x29 

108 

x22 

130 

x15 

87 

xl 

109 

x23 

131 

x16 

88 

x2 

110 

x24 

132 

x17 

TABLE  II   (Continued) 
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Element 

Variable 

Element 

Variable 

Element 

Variable 

133 

*18 

150 

"13 

167 

y12 

134 

x19 

151 

*21 

168 

y13 

135 

x20 

152 

y22 

169 

y21 

136 

x21 

153 

y23 

170 

y22 

137 

x22 

154 

y31 

171 

y23 

138 

x23 

155 

y32 

172 

y31 

139 

x24 

156 

y33 

173 

y32 

140 

x25 

157 

y41 

174 

y33 

141 

x26 

158 

y42 

175 

y41 

142 

x27 

159 

y43 

176 

y42 

143 

x28 

160 

y71 

177 

y43 

144 

x29 

161 

y72 

178 

y71 

145 

FDB. 

162 

y73 

179 

y72 

146 

FDBt 

163 

y81 

180 

y73 

147 

FDB3 

164 

y82 

181 

y81 

148 

yll 

165 

"83 

182 

y82 

149 

yi2 

165 

yll 

183 

y83 
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B.  Part  II 

Integration  of  the  equations  of  motion  requires  the  use  of  the  G  matrix 
from  Part  I,  and  that  the  initial  displacements  and  velocities  of  the  system 
are  known.  The  G  matrix  is  read  from  the  disk  on  which  it  was  stored.  The 
assumed  initial  conditions,  step  size,  and  number  of  steps  to  be  taken  are 
input  via  the  card  reader.  The  order  in  which  these  quantities  are  read  is 
shown  in  Table  III. 

The  last  ninety-seven  columns  of  G  are  used  in  a  new  matrix  H  which  takes 
on  the  form  of  the  right  hand  side  of  equation  (148) .  The  matrix  H  and  the 
initial  conditions  are  transferred  to  the  Runge-Kutta-Gill  subroutine.   This 
subroutine,  along  with  a  derivative  subroutine,  integrates  the  equations  for 
the  prescribed  number  of  steps.  The  resulting  displacements  and  velocities 
are  either  printed  out  or  stored  on  magnetic  disk  for  later  use.  Listings 
of  the  integration  routine,  Runge-Kutta-Gill  subroutine,  and  derivative 
subroutine  are  in  Appendix  A. 

The  y^'s  and  y^'s  are  defined  in  the  derivative  subroutine.  For  motion 
over  a  hard  flat  surface,  all  y's  and  y's  are  set  to  zero.   If,  however,  the 
front  wheels  see  a  sine  wave  in  the  Z  direction,  all  y's  and  y's  are  zero 
except  y.,,  y.,,  y-,  and  y.,.  The  expressions  for  these  variables  will  be 

^13  =  y23  =  a  3^n   "'     and     ^13  ™  ^23  "  a"  cos  "  (.151   a>b 

where  a  and  fj  are  some  constants. 

The  drawbar  loads  are  also  defined  in  the  derivative  subroutine.   If 
there  is  no  drawbar  load,  all  FDB's  are  set  to  zero.   If,  for  example,  the 
load  in  the  Y  direction  is  constant,  and  the  load  in  the  X  direction  varies 
sinusoidally,  the  drawbar  forces  would  be  input  as 
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TABLE  III 

ORDER  IN  WHICH  THE  INITIAL  CONDITIONS  ARE 
INPUT  TO  THE  COMPUTER 


Card Parameters 


1  X]^  x2,  x3,  X4,  x5,  x6>  x7,  xg 

xg,  x10,  xn,  x12,  x13,  x14,  x15,  x16 

x17'  x18'  x19'  x20'  x21*  x22»  x23*  x24 
x25'  x26*  x27'  x28*  x29*  1'  x2*  x3 
x4>  x5>  x6>  x7>  x8>  x9>  X10'  xll 
x12'  x13'  x14*  x15'  x16'  x17'  x18'  x19 

x20'  *21*  *22'  *23'  x24'  x25>  x26'  x27 

8  x28,  x2g 

9  DELT  (Time  increment) 

10  NSTEP*  (Number  of  steps  to  be  taken) 


*  This  value  Is  input  in  13  format,  while  all  other  quantities  are  F10.4 
format. 
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FDB!  =  a  sin  St     ,     FDB2  =  Y    and  FDB3  =  0  (152  a,b,c 

where  a,  6  and  Y  are  constants. 

The  values  of  the  FDB's,  y's  and  f's  are  arranged  in  a  vector  Y.  The 
order  of  these  functions  in  the  Y  vector  are  given  in  Table  IV.  Equations 
(151  a)  and  (151  b)  are  input  into  the  derivative  subroutine  as 

Y(6)  =  a  sin  St      ,     Y(9)  =  a  sin  St 

(153  a,b,c,d 
Y(2A)  -  aS  cos  St     ,     Y(27)  -  aS  cos  St 

Y  is  initialized  in  the  main  program,  eliminating  the  need  to  define  it 
in  the  derivative  subroutine  when  simulating  motion  over  a  hard  flat  surface 
with  no  drawbar  load. 
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TABLE  IV 
ORDER  OF  THE  FORCING  FUNCTIONS  IN  THE  VECTOR  Y 

Y(l)   =  FDBX  Y(14)  =  y42 

Y(2)   =  FDB2  Y(15)  =  y^ 

Y(3)   =  FDB3  Y(16)  =  <Jn 

*(*)   ■  7U  Y(17)  =  y72 

t<5)  =  y12  Y(18)  =  y?3 

Y(6)  =  y13  Y(19)  =  y81 

W)      =  yn  Y(20)  =  yg2 

Y(8)   =  y22  Y(21)  =  y83 

Y(9)  =  y23  Y(22)  -  yn 

Y(10)  =  y31  Y(23)  =  y12 

Ydl)  =  732  Y(24)  =  y13 

V(12)  =  y33  Y(25)  =  y21 

*G3)  =  y41  Y(26)  =  y22 


?(27) 

=  723 

Y(28) 

=  *31 

Y(29) 

"  732 

Y(30) 

-  733 

Y(31) 

"  *41 

Y(32) 

"  *42 

Y(33) 

"  *43 

Y(34) 

=  hi 

Y(35) 

=  *72 

Y(36) 

=  hi 

Y(37) 

"  *81 

Y(38) 

'  *82 

Y(39) 

=  *83 

VIII.   PRESENTATION  OF  RESULTS 

To  demonstrate  the  programs,  the  motion  of  a  tractor  rolling  over  a  half 
sinusoid  bump  with  its  left  rear  wheels  is  simulated.  The  parameters  for  a 
hypothetical  tractor  are  listed  in  Table  V.  These  parameters  are  input  into 
TRA  or  TRB  which  derive  the  equations  of  motion  for  this  tractor.  A  listing 
of  the  computer  generated  equations  of  motion  for  this  system  is  given  in 
Appendix  C.  The  equations  generated  by  TRA  are  identical  to  those  generated 
by  TRB;  therefore,  only  one  set  of  equations  is  listed  in  the  appendix. 

The  bump  is  input  via  the  Y  vector  in  the  derivative  subroutine.   Since 
only  the  left  rear  wheels  traverse  the  bump,  V33,  y^t   Yyj   and  y,,  are  the 
only  non-zero  elements  of  Y.  The  equation  of  a  half  sinusoid  bump  such  as 
that  shown  in  Figure  9  is 

y  =  h  sin  (TTvt/1)  (154 

where  h  is  the  height  of  the  bump,  1  is  the  length  of  the  bump,  and  v  is  the 
forward  velocity  of  the  tractor. 


1      -    vt 


FIGURE  9 

THE  BUMP  BEING  TRAVERSED  BY  THE  TRACTOR 
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TABLE  V 
LIST  OF  PARAMETERS  FOR  A  HYPOTHETICAL  TRACTOR 


Mass  (lb-sec/ft  and 

Inertia  (lb- 

-ft/sec^)  Values 

Ml  = 

1.83 

J1X 

■ 

0.57 

J1Y 

= 

1.14 

J1Z 

= 

0.57 

M2  = 

1.83 

J2X 

■ 

0.57 

J2Y 

■ 

1.14 

J2Z 

= 

0.57 

M3  = 

11.33 

J3X 

■ 

11.33 

J3Y 

■ 

22.66 

J3Z 

= 

11.33 

M4  = 

11.33 

J4X 

- 

11.33 

J4Y 

" 

22.66 

J4Z 

■ 

11.33 

M5  = 

128 

J5X 

= 

375 

J5Y 

■ 

900 

J5Z 

- 

1050 

Mg  = 

5.1 

J6X 

= 

9.0 

J6Y 

■ 

2.0 

J6Z 

= 

10.0 

M?  = 

11.33 

J7X 

= 

11.33 

J7Y 

■ 

22.66 

J7Z 

= 

11.33 

Mg  = 

11.33 

J8X 

■ 

11.33 

J8Y 

- 

22.66 

J8Z 

- 

11.33 

Dimensions  (ft) 

DU  = 

0.25 

D14 

- 

1.12 

D22  = 

0.25 

D24 

■ 

1.12 

D31  = 

0.0 

D32 

= 

0.33 

D33 

= 

0.33 

D34 

- 

2.00 

D41  = 

0.0 

D42 

= 

0.33 

D43 

= 

0.33 

D44 

= 

2.00 

D51  = 

3.67 

D53 

= 

1.55 

D55 

■ 

2.83 

D56 

■ 

2.83 

D57 

»  0.07 

D62  = 

0.0 

D63 

= 

0.0 

D64 

■ 

0.0 

D99 

= 

2.21 

D71  = 

0.0 

D73 

= 

0.33 

D74 

- 

2.00 

D81  = 

0.0 

D82 

- 

0.33 

D84 

= 

2.00 

L37  = 

1.50 

L53 

■ 

2.00 

L54 

= 

2.00 

L84 

= 

1.50 

X  = 

2.33 

Y 

= 

0.0 

Z 

" 

2.67 

DBj^ 

= 

2.50 

DB3 

=  0.67 

67 


TABLE  V  (Continued) 


Spring  Rates  (lb/ft) 


Ku  ■  16,000.0  K12  -  10,700.0  K13  -  22,600.0 

K21  =  16,000.0  K22  -  10,700.0  K23  -  22,600.0 

K31  =  18,000.0  K32  =  11,900.0  K33  =  20,500.0 

K41  =  18,000.0  K42  -  11,900.0  K43  =  20,500.0 

K?1  =  18,000.0  K?2  -  11,900.0  K?3  -  20,500.0 

K^  =  18,000.0  Kg2  =  11,900.0  Kg3  -  20,500.0 


Damping  Coefficients  (lb-sec/ft) 


Cu  -  88. 

,0 

C21  -  88, 

,0 

C31  -  134, 

,0 

C41  -  134, 

,0 

C?1  -  134, 

,0 

C81  -  134, 

,0 

C12 

= 

25 

,0 

C22 

■ 

23 

,0 

C32 

- 

32 

.0 

C42 

■ 

32 

,0 

C72 

- 

32 

.0 

C82 

■ 

32 

.0 

C13  ■  186, 

,0 

C23  =■  186, 

,0 

C33  -  248, 

0 

C43  =  248, 

,0 

C?3  -  248, 

,0 

Cg3  =  248, 

,0 

Material  Properties 


EI  -  1.33  x  106  JG  =■  1.06  x  106 
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Using  equation  (153) ,  the  expressions  for  y,,  and  y7,  are  then 

Y(12)  =  Y(18)  =  y33  =  y?3  =  -.416  sin  (4.61t)  (155 

and  y,,  and  y,,  are  given  by 

Y(30)  -  Y(36)  =  y33  =  y?3  =  -1.914  cos  (4.61t)  (156 

for  the  interval  0  <_  t  £_  .6815,  and  are  zero  otherwise.  These  values  are 
negative  sinze  Z  is  positive  downward.  These  are  the  expressions  for  Y  used 
in  the  listing  of  the  derivative  subroutine  in  Appendix  A. 

The  assumed  initial  conditions  are  such  that  the  tractor  is  initially  at 
rest,  that  is  all  velocities  and  displacements  are  zero.  A  step  size  of  .001 
sec.  is  used,  and  1,000  steps  are  taken.  The  results  of  every  twentieth  step 
are  listed  in  Appendix  C. 

To  compare  the  results  of  this  model  with  those  of  a  model  having  rigid 
axles,  the  programs  are  run  again  using  the  parameter  listed  in  Table  V, 
except  for  the  values  of  EI  and  JG,  which  are  increased  by  a  factor  of  ten. 
Figures  10  and  11  show  the  verticle  deflections  of  the  centers  of  gravity  of 
the  left  rear  wheels,  for  both  the  soft  and  stiff  cases. 
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IX.   DISCUSSION  AND  CONCLUSIONS 

As  previously  mentioned,  the  equations  generated  by  TEA  and  TRB  are 
identical,  which  implies  that  they  are  in  fact  the  correct  equations  of 
motion  for  the  system  provided  the  original  assumptions  are  valid.  Both 
of  these  programs  have  certain  advantages  as  well  as  drawbacks. 

TRA  is  relatively  expensive  to  use,  since  it  involves  the  reduction  of 
an  eighty-six  by  one  hundred  eighty-three  matrix.  To  compile  this  program, 
execute  it,  and  print  the  equations  of  motion  on  both  paper  and  magnetic 
disk  requires  approximately  fifty-five  seconds  of  CPU  time  on  the  Itel  AS/5 
computer.  The  first  fifty-seven  rows  of  the  reduced  E  matrix,  in  this  pro- 
gram, are  expressions  for  the  forces  and  moments  acting  on  the  system  in 
terms  of  x's,  x's  and  the  forcing  functions.  These  expressions  may  be  used 
to  evaluate  the  force  acting  on  an  axle,  for  example,  at  some  instant  in 
time. 

When  simulating  a  tractor  with  very  stiff  axles,  TRA  requires  some 
equations  to  be  normalized  with  respect  to  a  very  small  quantity.  This 
introduces  some  round  off  error  which  may  be  multiplied  through  the  reduc- 
tion process.  The  result  is  a  few  unwanted  variables  being  multiplied  by 
quantities  of  the  order  of  101  or  less,  while  the  correct  variables  are 
being  multiplied  by  coefficients  of  order  10  or  better.  These  unwanted 
variables  have  negligable  effect  on  the  integration  of  the  equations  of 
motion. 
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TRB  requires  the  reduction  of  a  twenty-nine  by  one  hundred  twenty-six 
matrix,  which  makes  it  much  faster  to  execute  than  TRA.  To  compile,  execute, 
and  print  the  equations  of  motion  with  this  program  requires  about  twenty- 
three  seconds  of  CPU  time,  less  than  half  that  of  TRA. 

This  program  does  not  yield  expressions  for  forces  or  moments  acting 
on  the  system,  but  if  these  expressions  are  not  needed,  the  savings  in  CPU 
time  is  substantial.  Round  off  and  inherited  error  are  also  reduced  with 
this  program,  since  it  requires  fewer  mathematical  operations. 

The  CPU  time  required  to  execute  the  integration  routine  depends  on  the 
number  of  steps  taken,  which  depends  on  the  time  increment  used.  The  size 
of  the  time  increment  is  not  completely  arbitrary.  Too  large  a  step  size 
will  cause  the  system  of  equations  to  oscillate  wildly  or  blow  up,  and  too 
small  a  step  size  requires  that  excessively  many  steps  be  taken.  The  proper 
step  size  must  be  determined  by  trial  and  error.  As  a  rule  of  thumb,  stiffer 
axles  dictate  a  smaller  step  size.   CPU  time  may  also  be  saved  by  not  print- 
ing values  after  every  step. 

For  a  given  system,  the  amount  of  time  needed  to  perform  a  step  does 
not  depend  on  the  step  size  to  any  noticeable  degree.  For  the  example  of 
section  VIII,  one  thousand  steps  with  a  step  size  of  .001  sec,  printing 
every  tenth  step,  required  approximately  twenty-five  seconds  of  computer 
time. 

In  deriving  the  equations  of  motion  by  the  energy  method,  it  should  be 
noted  that  the  potential  energy  in  the  deflected  axles  is  derived  by  the 
method  given  In  Appendix  B.  Most  of  the  previous  models  have  had  rigid 
axles  which  eliminate  this  energy  term  altogether.  Mather's  model  had 
elastic  axles,  but  his  expression  for  the  potential  energy  of  a  deflected 
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section  of  the  axle  was 

U  =  h   (3EI/L3)  xt2  (157 

This  implies  no  moment  is  applied  to  the  rear  axle,  which  is  not  really  a 
valid  assumption.  Use  of  equation  (157)  produced  a  set  of  equations  that 
were  radically  different  from  those  generated  using  equation  (B  12) . 

Figures  10  and  11  show  that  there  is  noticeable  difference  between  this 
model  and  models  that  treat  the  rear  axles  as  rigid  members.  Another  simu- 
lation was  attempted  in  which  the  parameters  of  Table  V  were  used,  except 
that  EI  and  JG  were  multiplied  by  one  hundred.  The  resulting  system  was  so 
stiff  that  the  step  size  needed  to  successfully  integrate  the  equations  of 
motion  was  prohibitively  small. 

By  allowing  the  rear  axles  to  deflect,  the  stress  in  an  axle  for  a  given 
deflection  could  be  determined.  The  knowledge  of  this  stress  level  would 
permit  one  to  predict  when  the  axle  would  fail. 

In  conclusion,  the  computer  programs  presented  here  are  a  very  general 
base,  with  which  most  any  dual  rear  wheeled  tractor  can  be  simulated.  The 
user  is  allowed  a  good  deal  of  freedom  in  selecting  tractor  parameters,  and 
an  unlimited  variety  of  forcing  functions  via  the  subroutine  DERIV.  Param- 
eters and  initial  conditions  are  easily  varied  to  make  comparisons  between 
different  systems  or  different  situations. 

This  work  leads  to  a  vast  quantity  of  further  research.  Among  this 
research  is  experimental  verification  of  this  model,  extension  of  the 
computer  program  to  determine  stress  levels  in  the  rear  axles,  and  con- 
sideration of  the  tires  as  nonlinear  springs. 
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APPENDIX  A 


LISTINGS  OF  THE  COMPUTER  PROGRAMS 
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THIS    PROGRAM    GENERATES    THE    EQUATIONS    OF    MCTICN    USING   THE    EOUAT IONS 
DEVELOPED    IN    SECTION    5A. 

IMPLICIT  REAL(A-Z) 
INTEGER  IfJ,K.N,Nl 
DIMENSION    A(o6,29) 

DIMENSION    8(66,29) 
DIMENSION    C186.29) 
DIMENSION    D(86,39) 
DIMENSION    E186.183) 
DIMENSION    FI86.57)      ' 
DIMENSION   G(29,126) 
DIMENSION    NAMEU26) 
REWIND    10 

THE    NEXT    15    CARDS    READ    IN   THE    VARIOUS    TRACTOR    PARAMETERS. 


TRA 

1 

TRA 

2 

TRA 

3 

TRA 

4 

TRA 

5 

TRA 

6 

TRA 

7 

TRA 

8 

TRA 

9 

TRA 

10 

TRA 

11 

READ  (5,180)  Ml.JIX.JlY, J1Z,D13,D14 
READ  (5,1801  M2.J2X.J2Y,  J2Z, 022,024 
REAO    (5,180)    M3.J3X,  J3V,  JSZ.C^' 


READ    (5,180)    L37,L53,L54,L84 
READ    (5,200)     EI.GJ  TRA 

READ    (5,180)     XSAR,YBAR,Z6AR,DB1,D83  TRA 

READ    (5,180)    KU,K12,K13,K21,K.22,K23,K31,K32,K33,K41,K42,K43,K71,K.TRA 

*i2,K73,K61,K82,KE3  TRA 

READ    (5,180)    Cli,C12,C13,C21,C22,C23,C31,C32,C33,C41,C42,C43,C71,CTRA 

*72,C73 ,C81,C32,C83  TRA 


C  N    IS    ThE    NUMBER    OF    ROWS    AND    Nl    THE    NUMBER    OF    COLUMNS    OF    THE    MATRIX    E. 

C 

N   =    86 

M    =    183 
C 

C      THE  NEXT  11  CARDS  INITIALIZE  THE  A,8,C,0  AND  F  MATRIC1ES. 
C 

DO  40  1=1,86 

00  10  J=l,29 

AU.J)  «  0.0 

8(1,  J)  ■  0.0 


TRA 

38 

TRA 

39 

TRA 

40 

TRA 

41 

TRA 

42 

:  E. 

TRA 

43 

TRA 

44 

TRA 

45 

TRA 

46 

TRA 

47 

TRA 

48 
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C(I.J)  =  0.0 

10 

CONTINUE 

DO  20  J=l 

57 

20 

F(I,J)  ■  0.0 

DO  30  J=l 

39 

30 

ci  i, j)  ■  o.o 

40 

c 
c 

f 

CONTINUE 

THE  NEXT  79  CARDS  EVALUATE  T 

u 

AI39.1)  ■ 

Ml 

A139.5)  = 

M1*(D53+D62) 

A(39,6)  = 

Ml*(D99+D13+YaARJ 

AI40.2)  = 

Ml 

A140.4)  ■ 

-M1*D53 

A(40,6)  = 

Ml»051 

A (40,7)  ■ 

-Ml*062 

A141.3)  * 

Ml 

A(41,4>  => 

-M1*YBAR 

A(41,5)  = 

-M 1*051 

A141.7)  = 

-M1*(C99+013) 

AI42.7)  = 

J IX 

A143.8)  = 

J1Y 

A<44,6)  = 

J1Z 

A(4S,1)  = 

M2 

A(45,5)  = 

M2*(D53+C62)  . 

A(45,6)  = 

-M2*(D99+022-YeAR) 

A(46,2)  = 

M2 

A(46,4)  = 

-M2*053 

A(46,6)  = 

M2*D51 

A (46, 7)  = 

-M2*D62 

A(47,Jl  ■ 

M2 

A<47,4)  = 

-M2*YBAR 

A147.5)  = 

-M2+051 

A(47,7)  = 

M2*(O99+022) 

A(48,7I  = 

J2X 

A(49,9)  = 

J2Y 

A(50,6)  = 

J2Z 

A(51,14) 

=  M3 

A(52,2)  = 

M3 

A(52,4)  = 

-M3*D57 

A(52,6)  = 

-M3*XEAR 

A(53,10) 

=  M3 

AI54.26) 

=  J3X 

A<55,13) 

=  J3Y 

AI56.22) 

=  J3Z 

A157.15I 

=  M4 

A( 5q,2J  = 

M4 

A(58,4)    =    -M4*057 
A(58,6)    =    -M4*X8AR 
A( 59,11)    =    M4 
A160.27)    =    J4X 
AU1.19)    ■    J4Y 
AI62.23)    •    J4Z 
A (63,1)    =    M5 
A(64,2)    =    M5 
A(c5,3)    =    M5 
A166.4)    =    J5X 
A167.5)    =    J5Y 
A(63,6)    ■    J5Z 


TRA 

49 

TRA 

50 

TRA 

51 

TRA 

52 

TR,' 

53 

TRA 

54 

TRA 

55 

TRA 

56 

TRA 

57 

TRA 

sa 

TRA 

59 

TRA 

60 

TRA 

61 

TRA 

62 

TRA 

63 

TRA 

64 

TRA 

65 

TRA 

66 

TRA 

67 

TRA 

68 

TRA 

69 

TRA 

70 

TRA 

71 

TRA 

72 

TRA 

73 

TRA 

74 

TRA 

75 

TRA 

76 

TRA 

77 

TRA 

78 

TRA 

79 

TRA 

80 

TRA 

81 

TRA 

82 

TRA 

83 

TRA 

84 

TRA 

85 

TRA 

86 

TRA 

87 

TRA 

88 

TRA 

89 

TRA 

90 

TRA 

91 

TRA 

92 

TRA 

93 

TRA 

94 

TRA 

95 

TRA 

96 

TRA 

97 

TRA 

98 

TRA 

99 

TRA 

100 

TRA 

101 

TRA 

102 

TRA 

103 

TRA 

104 

TRA 

105 

A(69 

1) 

• 

M6 

A(69 

5) 

= 

M6*( 053+063) 

A(69 

6) 

• 

M6*YEAR 

A(70 

2) 

» 

M6 

A170 

4) 

m 

-M6*D53 

A(70 

61 

s 

Mo* 051 

A(70 

7) 

= 

-M6*063 

A(71 

3) 

■ 

M6 

79 


TRA  106 

TRA  107 

TRA  108 

TRA  109 

TRA  110 

TRA  111 

TRA  112 

TRA  113 

A(71,4)  ■  -M6*YBAR  TRA  114 

A(71,5)  =  -M6*051  TRA  115 

A(72,7)  =  J6X  TRA  116 

AI73.5)  =  J6Y  TRA  117 

A174.6)  =  J6Z  TRA  118 

A( 75,16)  =  M7  TRA  119 

AI76.2)  =  M7  TRA  120 

A176.4)  =  -M7*057  ,                        TRA  121 

A(7b,6)  =  -H7*XBAR  TRA  122 

AI77.12)  ■  M7  TRA  123 

A(78,28)  =  J7X  TRA  124 

A( 79,20)  =  J7Y  TRA  125 

A(80, 24)  ■  J7Z  TRA  126 

A(Sl.U)  =  M8  TRA  127 

A(82,2)  =  M8  TRA  128 

A(82,4)  =  -M8*057  TRA  129 

A(82,6)  =  -M8*XBAR  TRA  130 

A(83,13)  =  M8  TRA  131 

A184.29)  =  J8X  TRA  132 

A( 85,21)  =  J8Y  TRA  133 

A186.25)  =  J8Z  TRA  134 
C 

C      THE  NEXT  52  CARDS  EVALUATE  THE  NCNZERO  ELEMENTS  CF  THE  MATRIX  e. 
C 

8(21,1)  =  Cll  TRA  135 

8(21,51  =  Cll*(D53*Dfc2)  TRA  136 

6(21,6)  =  C11*(C99+D13-»YEAR)  TRA  137 

5(21, d)  =  (.11*014  TRA  138 

8(22,2)  =  C12  TRA  139 

3(22,4)  =  -C12*053  TRA  140 

B(22,6)  =  C12*051  TRA  141 

E122.7)  -    -C12*(C62+C14)  TRA  142 

B(23,3)  =  C13  TRA  143 

B(23,4)  =  -C13*YBAR  TRA  144 

8(23,5)  =■  -C13*D51  TRA  145 

B(23,7)  =  -C13*(D99+C13)  TRA  146 

BI24.1)  =  C21  TRA  147 

BI24.5)  =  C21*(C53+Q62)  TRA  148 

B(24,6)  =  -C21*(D99+022-YBAR)  TRA  149 

£124,9)  =  C21*024  TRA  150 

8(25,2)  =  C22  TRA  151 

3(25,4)  =  -C22*D53  TRA  152; 

E(25,6)  =  C22*051  TRA  153 

BI25.7)  =  -C22*(D62+D24)  TRA  154 

6(26,3)  *  C23  TRA  155 

B(2o,4)  =    -C23*Y8Afl  TRA  156 

B126.5)  =  -C23*D51  TRA  157 

3(26,7)  =  C23*(C99+022)  TRA  158 

BI27.14)  =  C31  TRA  159 

e(27,18)  =  C31*D34  TRA  160 

E(28,2)  =  C32  TRA  161 

BI26.4)  =  -C32*057  TRA  162 
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• 

£(28,6)    =    -C32*XBAR 

B<28,26)    =   -C32*034 

B(29,10)    =   C33 

8(30,15)    =   C41  ' 

BI30.19)    =   011*044 

8131,2)    =    C42 

B(31,4)    =    -C42*D57 

BI31.6)    =    -C42»XBAR 

6(31,27)    =    -C42*044  ! 

B(32,ll)    =  C43 

B(33,16)    ■   C71 

E(33,20)    =   C71*074 

BI34.2)    =   C72 

B134.4)    =   -C72*C57 

E(34,6)    -   -C72*XBAR 

B(34,2£)    =    -C72*D74  '  \ 

B(35,12)    =    C73 

B(  3o,17)    =   C31 

6(36,21)    =    CS1»084 

BI37.2)    =    C82 
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F( 65, 5<t) 

» 

-1 

Ff  66,  281 

■ 

-(L64+082) 

F(86,57) 

= 

-081 

F( 66,55) 
C 

C              THE    NEXT 
C 

00    80    1=1 

* 

-1 

iC 

1    CARDS    CAS1 

,86 

DO    50    J  =  l 

,57 

50    Ed.Jl    = 

F(I,J) 

DO    60    J=l 

,29 

E(I ,J+57) 

:    -Ad, J) 

El  I.J  +  66) 

'   CII.JI 

60    E(I,J+115) 

■    8(1, J) 

DO    70    J=l 

,39 

70    ElI.J+144) 

=    -0(1, JJ 

80    CONTINUE 

CALL    GAUSS 

(E,N,NU 

86 

TRA    514 

TRA    515 

TRA    516 

TRA    517 

TRA    518 

TRA    519 

TRA    520 

TRA    521 

TRA    522 

TRA    523 

TRA    524 

TRA    525 

TRA    526 

TRA    527 

TRA   528 

TRA    529 

TRA    530 

TRA    531 

TRA    532 

TRA    533 

TRA    534 

TRA    535 

TRA    536 

TRA   537 

TRA    538 

TRA    539 

TRA    540 

TRA    541 

TRA    542 

TRA    543 

TRA   544 

TRA    545 

TRA    546 

TRA    547 

TRA    548 

TRA    549 

TRA    550 

TRA    551 

TRA  552 
TRA  553 
TRA  554 
TRA  555 
TRA  556 
TRA    557 


A,B,C,0    AND    F    MATR1CIES    INTO    THE    MATRIX 


C 

C  THE    NEXT    4    CAROS    CAST    THE    LAST    126    COLUMNS    OF    THE   LAST    29   ROWS    OF 


(    E. 

TRA 

558 

TRA 

559 

TRA 

560 

TRA 

561 

TRA 

562 

TRA 

563 

TRA 

564 

TRA 

565 

TRA 

566 

TRA 

567 

TRA 

568 
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C      ThE  MATRIX  E  INTO  ThE  MATRIX  G. 
C 

DO  100  1=1,29  TRA  569 

DO  SO  J=l,126  TRA  570 

90  G(I,J)  ■  E(I  +  57,J<-57)  TRA  571 

100  CONTINUE  TRA  572 

WRITE  (10)  G  TRA  573 

C 

C  THE   NEXT    13    CARDS    PRODUCE   A    FORMATTED   COPY    OF    THE    NCNZERO    ELEMENTS 

C  OF    THE    MATRIX    G. 

C 

DO    120    1=1,29  TRA  574 

WRITE    (6,190)     I  TRA  575 

DO    110   J=l,126  TRA  576 

IF    (J.Efl.l)    ID   =    1  TRA  577 

IF    (AES(G(I, J)) .LE. 0.0001)    GO    TO    110  TRA  578 

IF    (IC.EQ.l)    WRITE    (6,130)    G  ( I  ,  J)  ,  NAPE  (  J)  TRA  579 

IF    (iC.Ed.2J    WRITE    (6,140)    G ( I , J ), NAME ( J )  TRA  580 

IF    (10. Eg. 3)     WRITE    (6,150)    G ( I , J ) , NAKE ( J)  TRA  581 

IF    (IC.EC.4)    IrlRITE    (6,160)    G(I,J),NAMEU)  TRA  582 

IF    (10.EQ.4)    ID   =    1  TRA  583 

Id   =    ID+1  TRA  584 

110    CONTINUE  TRA  585 

120    CONTINUE  TRA  586 

STOP  TRA  587 

130  FORMAT  (•  •,20X,F3.1, IX, A4,'=')  TRA  588 

140  FORMAT  ( '+',T31,'  +  ('  ,F13.4,'  )'  ,A4)  TRA  589 

150  FORMAT  (  '  +  '  ,T51,  '+  (  '  ,F  13  .4, '  )  '  ,A4)  TRA  590 

160  FORMAT  ( ' + ' ,T71 , ' + ( ' ,F13 .4 , • ) • , A4, / )  TRA  591 

170  FORMAT  (20A4)  TRA  592 

180  FORMAT  (8F10.2)  TRA  593 

1*0  FORMAT  (  '-'  ,47X,'ECUATI0M  ,15)  TRA  594 

200  FORMAT  (2E10.4J  TRA  595 

END  TRA  5S6 
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THIS    PROGRAM    GENERATES    ThE    EQUATIONS   OF   MOTION    USING    THE    EQUATIONS 

C 

IMPLICIT    REAL(A-Z) 

INTEGER    I,J,K,N,N1 

GIMENSICN    AI29.29) 

DIMENSION    B(29,29) 

DIMENSION    0(29,291 

DIMENSION    0129,39) 

DIMENSION    G(29,126) 

DIMENSION    NAME! 126) 

REWIND    10 
C 
C  THE   NEXT    15    CARDS    READ    IN    THE    VARIOUS    TRACTOR    PARAMETERS. 


C 


c 


TRS 

1 

TRB 

2 

TRB 

3 

TRB 

4 

TRB 

5 

TRB 

6 

TRB 

7 

TRB 

3 

TRB 

9 

READ    (5,140)  M1,J1X,J1Y, J1Z, 013,014  TRB  10 

READ    (5, 140)  M2.J2X, J2Y, J2Z,022,024  TRB  xl 

READ    (5,140)  M3.J3X.J3Y,  J3Z,  031,032, 033,034  TRB  12 

READ    (5,140)  M4.J4X, J4Y, J4Z, 04 1 ,042 ,043 ,044  TRB  13 

REAO    (5,140)  M5,J5X,J5Y, J5Z, D51, 053 , 055, D56, 057  TRB 


WRITE     (6,140)     C11,C12,C13,C21,C22,C23,C31,C32,C33,C41,C42,C43,  C71.TRB 

*C72,C73,Cfil,C82,Ce3  TRa 

READ  (5,130)  INAMEU),  1=1,126)  TRB 

N  IS  THE  NUMBER  OF  ROWS  ANO  Nl  THE  NUMBER  CF  COLUMNS  OF  THE  MATRIX  G. 

N  =  29 

Nl  =  126 
C 
C      THE  NEXT  10  CARDS  INITIALIZE  THE  A.B.C  ^N0  0  MATRICIES. 


A(I,J)    =    0.0 
e(I,J)    =0.0 

ci i , j)  »  o.o 

10    CONTINUE 


14 


READ  (5,140)  N6.J6X, J6Y, J6Z, 062,063, 064,099  TRB  15 

READ  (5,140)  M7.J7X.J7Y, J7Z, 071,073, 074  TRB  16 

READ  (5,140)  MS.JdX,  JOY,  J6Z,0ai,0o2,C84  TRB  17 

READ  (5,140)  L37,LS3,L54,L84  TRB  la 

REAO  (5,160)  EI.GJ  TR8  ll 

READ  (5,140)  XBAR,YBAR.,ZEAR,CB1,DB3  TRB 


REAO    (5,140)    KU,K12,Ki3,K21,K22,K23,K31,K32,K33,K41,K42,K43,K71,KTRB  21 

*72,K7S,K81,K82,K83                                                                                                                           TRB  22 

READ    (5,140)    C11,C12,C13,C21,C22,C23,C31,C32,C33,C41,C42,C43,C71,CTRB  23 

*72,C73,C81,C82,C83                                                                                                                           TRB  24 

WRITE     (6,140)    M1,J1X,J1Y,J1Z,013,014                                                                             TRB  25 

WRITE    (6,140)    M2,J2X,J2Y,J2Z,022,024                                                                             TRB  26 

WRITE     (6,140)    M3,J3X,J3Y,J3Z,  031, 032,033,034                                                         TRB  27 

WRITE     (6,140)     M4,J4X,J4Y,J4Z,041,042,D43,C44                                                             TRB  28 

WRITE    (6,140)    M5,J5X,J5Y,J5Z, 051, 053,055, 056, 057                                               TRB  29 

WRITE     (6,140)     M6.J6X.J6Y ,J6Z, J62 , J63 ,064, 099                                                              TRB  'O 

WRITE    (6,140)    M7,J7X,J7Y,J7Z,071,D73,074                                                                   TRB  31 

WRITE     (6,140)    Ma,J8X,J8Y,J8^,D81,D82,0a4                                                                         TRB  32 

»RITE    (6,140)    L37,L53,L54,L84                                                                                               TRB  ** 

WRITE    (6,160)    EI.GJ                                                                                                                         TRB  34 

WRITE    (6,140)    XEAR, YEAR, ZBAR, 081, 0B3                                                                             TR8  35 

«RITE     (6,140)     K11,K12,K13,K21,K22,K23,K31,K32,K33,K41,K42,K43,K71,TRB  36 

*K72,K73,K81,K82,K83                                                                                                                        TRB  37 

38 
39 
40 


TRB       41 
TRB      42 


CO    30     1=1,29  TRfl       ., 

DO    10    J=l,29 


TRB  44 

TRB  45 

TRB  46 

TRB  47 

TRB  48 
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00    20    J=l,39                                                                                                                                          TDO  .. 

D(I.J)    .0.0                                                                                                                                          }RB  49 

20    CONTINUE                                                                                                                                                     '„  *° 

30    CONTINUE                                                                                                                                                     ™J  *} 

£                                                                                                                                                                                                                                         1 Kd  3  2 

C  THE    NEXT    16   CARDS    EVALUATE    CONSTANTS   TO    BE   USED    LATER. 

KA1    =     <12.*EI)/L53**3.                                                                                                                 T0B  -, 

KA2    =    (12.*EI)/L37**3,                                                                                                                 „,  If 

KA3   =    U2.*EI)/L54**3.                                                  ■                                                                {.*"  |* 

TRB  56 

TR8  57 

TRB  58 


KA4    =     <12.*EI)/L84**3. 

KA5   =    <6.*EI)/L53**2. 

KA6    =    16.*EI)/L37**2. 

KA7    =     (6.*£I)/L54**2. 

KA8    «     <6.*EI)/L84**2. 

KA9    •     14.*EI)/L53 

KA10   =    (4.*EI)/L37 

KA11    =    (4.*EI)/L54 

KA12    =    (4.*EI)/L84 

KA13    -    GJ/L53  ™8      64 


TRB  59 

TR8  60 

TRB  61 

TRB  62 

TRB  63 


TRB       66 
TRB      67 


KA14    ■    GJ/L37 
KA15    =    GJ/L54 

c  KA16    "    GJ/L8*  TRB      63 

C  THE    NEXT    58   CAfiOS    EVALUATE    THE    NONZERO    ELEMENTS    CF    THE    MATRIX    A. 


TRB  69 
70 
71 
72 
73 
74 
75 

TRB  76 
77 

TRe      78 


A(l,l)    =    M1+M2+M5+M6 

All, 5)    =    D53*(M1+M2+M6)+062*(M1+M2I+C63*M6  trb 

AI01.06J    =    YbAR*<Ml-M2+M6)+099*(Ml-M2)+D13*Ml-022*M2  TRB 

A(02,02)     =    M1+M2+M3+M4+M5+M6+M7+M6  TRB 

A(02,04)    ■    -053*<M1+M2*-M6)-057*(M3  +  C4+M7+H8)  TRB 

A102.06)    «    051*lMltM2+M6)-XBAR*(M3+M4+M7+M8j  T»B 

AIC2.C7)     =    -062*(^l  +  M2)-0o3*M6                                                                                                    ,„  JT 
AI03.03)    =    M1+M2+M5+M6                                                                                                                          -Irr 
AI03.04)    =    -YEAR*(M1+M2+M6)                                                                                                             Jon 
A(03,05)    =    -D51*(M1  +  M2+M6I                                                                                                                jj! 

AI03.07)    =    D99*(M2-M1)-D13*M1  +  D22*M2                                                                              too  7„ 

AI04.02)     =    -D53»IMl+M2+Mo)-D57»IM3+C4+M7+«8J                                                              TRB  in 

AI04.03)    =    -YBAR»(M1+M2+M6)                                                                                                             J™  8° 

A(04,04i    =     ID53**2-.YEAR«*2)*IM1+M2+M6)+057**2*M3+M4+M7«-MS)  +  J5X         TR3  82 

AI04.C3)     =    YBAR*D51*(M1+M2+M6J                                                                                                     iZZ  f* 

AI04.06)    =    -051*Q53*(M1+M2+M6)+XBAR*D57*(M3+M4+M7+M8)                                      TRB  84 

/!B!a'2nL;„?53*062*(M1+M2J+i;53*06j*M6+YaAR*099*(Mi-M2'  +  'BAR*013*MlTRa  65 

»-rt)AK*U22*M2                                                                                                                                                      T0B  - 

AI05.01)    =    lD53*Cb2i*lHl*HU*lOS3*atii*H6                                                            trb  at 

A(05,03)     =    -D5HMM1+M2+M61                                                                                                                  „„  SI 

A(05,C4)    =    D51*YBAR*(M1+M2+M6)                                                                                                     j„  f| 

*A,i1i',1V     =    ID53  +  062)**^*<M1+M^1+I05^D63)**2»M6+D51»*2*(M1+M2*M6)*TRB  90 

Ton  oi 

*3j2mi»A6    'D"  +  062)*(lYBAR+D99+013,*Mlt('aAR-0«9-022J*M2)*(D53+D6TRB  92 

A(05,C7)    =    D5i*((D99+D13)*Ml-(099+D22)*M2)                                                                   TRB  94 

A(06,01)    .     <YBAR+D99+O13J*MltlYBAR-099-022)*M2+YflAR*M6                                  TRB  95 

AIC6,0^1     =    C51*(Ml+M2*-M6)-XBAR«(f<3  +  ><4+M7  +  Me)                                                              TRB  96 

a!^'?=!    =    7D5i*D53*(M1+M2'-051*053*M6+057»XBAR*lM3+M4-M7+MSI               TRB  97 

•  MM^»!.J2  53  +  C62""^ARtC99>£13)*Ml+(YBAR-D99-D22)*M2)  +  YBAR*(0TRB  98 

A106.C6)  =  (YBAR  +  D99)  ™"  " 
A106.C6)     =    <Y8Afi  +  D99  +  D13)*»2*Ml  +  lYEAR-C99-C221**2*M2+D51**2*IMl+M2TR8    101 

*+M6)+XEAR».2.<M3+M4+M7-M8)+YBAR«2*M6*JlZ  +  J2Z+J5ZtJ6Z  TRB    102 
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A(06,07) 

■ 

-051*L62*(M1+M2)-051*C63*M6 

TR8    103 

A(C7,C2) 

= 

-D62*(M1+M2)-063*M6 

TRB    104 

A107.03) 

■ 

(D99+013)*(M2-M1) 

TRB    105 

A  ( 07 1 0*1 

= 

D53»062*(MltM2)+063*05j*M6+YBAR»(C99+013)*(Ml-M2) 

TRB    106 

A(07,05) 

■ 

( 099*013) *0  51»(M1-M2) 

TRB    107 

A<07,06) 

* 

-062*D5  1*(M1+M2) -063*05 1*M6 

TRB    108 

A(07,C7) 

= 

062**2*  (Ml+M21  +  063**2*M6+(  099  +  022  )**2*(M1+M2)+J  IX +J2X+JTRB    109 

*6X 

TRB    110 

A(08,C8) 

= 

J1Y 

TRB    111 

A(09,09) 

= 

J2Y 

TRB    112 

A(10, 10) 

= 

M3 

TRB    113 

A(  11,11) 

= 

M4 

TRB    114 

AI12.12) 

■ 

M7 

TRB    115 

A(  13 ,13) 

■ 

M8 

TRB    116 

A(  14,14) 

= 

H3 

TRB    117 

A(1S, 15) 

> 

M4 

TRB    118 

A116.16) 

= 

M7 

TRB    119 

A(17,17) 

■ 

M3 

TRB    120 

AU8.18) 

■ 

J3Y 

TRB    121 

A(  19,19) 

■ 

J4Y 

TRB    122 

A(20,20) 

■ 

J7Y 

TRB    123 

A<21,21) 

= 

jaY            ; 

TRB    124 

AI22.22) 

= 

J3Z 

TRB    125 

AI23.23) 

= 

J4Z 

TRB    126 

A(24,24) 

= 

J7Z 

TRB    127 

A(  25,25) 

= 

J8Z 

TRB    128 

A(26,26) 

a 

J3X 

TRB    129 

A(27,27l 

= 

J4X 

TRB    130 

A(28,28) 

= 

J7X 

TRB    131 

c 

c 
c 

A129.29) 

■ 

J8X 

TRB    132 

THE   NEXT 

97   CARDS    EVALUATE    THE    NGNZERO    ELEMENTS    CF    THE    MATRIX    E. 

BlOl.Cl) 

= 

C11+C21 

TRE    133 

6101,05) 

= 

(C53+D62)*(Cll+C21) 

TRB    134 

B101,C6) 

■ 

C11*(Y6AR+099+013)+C21*(Y6AR-099-C22) 

TRB    135 

e(01,C8) 

= 

011*014 

TRB    136 

6(01,09) 

* 

C21*C24 

TRB    137 

B< 02,02) 

m 

C12+C22+C32+C42+C72+C82 

TRB    138 

EI02.04) 

» 

-D53*(C12+C22)-057*(C32+C42+C7  2+C82) 

TRB    139 

6(02,06) 

= 

051*(C12+C2  2)-X6AR*(C32+C42+C7  2+C82) 

TRB    140 

B(02,C7) 

s= 

-C 12* I 062+J 14 )-C22*( 062+0241 

TR8    141 

6(02,26) 

■ 

-C32»C34 

TRB    142 

8(02,27) 

= 

-042*044 

TRB    143 

6<  02>28) 

= 

-072*074 

TRB    144 

6(02,29) 

= 

-C82*C84 

TRB    145 

BI03.03) 

= 

C13+C23 

TRB    146 

6(03,04) 

= 

-YEAR*(C13+C22) 

TRB    147 

6(03,05) 

■ 

-D51*(C13+C23) 

TRB    148 

8(03, C7) 

■ 

-C 13* (0  99+0 13 )+C23*( 099*022) 

TR8    149 

6(04,02) 

= 

-053*(C12+C22)-05  7*(C32+C42+C7  2+Ce2) 

TRB    150 

6(04,03) 

= 

-YEAR*(C13+C23) 

TRB    151 

8(04,04) 

= 

T53**2*(C  12+C22)+YBAR**2*(C13+C23)+057**2*(C32+C42+C72 

+TRB    152 

*C82) 

TRB    153 

6( 04,05) 

■ 

D51*YBAR*(C13+C23) 

TRB    154 

6(04,06) 

m 

-C51*05  3*(C12+C22)+XbAR*0  57*(C22+C4  2+C72+C82) 

TRB    155 

81  04,07) 

= 

053* ( (062+014)*C12+(D62+C24)*C22)+Y6AR*( (099  +  013 )  *C  13- 

(TRB    156 

*C99+D22)*C23) 

TRB    157 

E(04,26) 

= 

032*057*034 

TRB    158 

8(04,27) 

■ 

042*057*044 

TRB    159 
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£(04,28)  =    C72*C57*D74  TRB    160 

3(04,29)  =   C32*D57*DS4  TR8    161 

e(05,01)  =    (D53  +  D62)*(C11+C21)  TR8    162 

B(05,03)  =   -051*(C13+C23)  TR8    163 

B(C5,04)  ■    YBAR*D51*(C13+C23)  TRB    164 

E(05,05)  =    (D53+062)**2*(i.ll+C21)»D51**2»lC13+C23)  TRB    165 

BIOS, Co)  ■    (D53+062)*( (Y8AR+G99+D13) *C 11+ ( YBAR-D99-C22)*C21 )  TRB    166 

8(05, C7)  =   051*( (099+0 13)*C13-(099+D22)*C23)  TRB    167 

fi(05,08)  ■    Cll»C14*(£53+C62)  TRB    168 

B(C5,'C9)  =  C21*C24*(053+£62)  TRB    169 

B(06,01)  =    Cll*(YbAR+099+D13)+C21*(YBAR-D99-022)  TRB    170 

B(06,02)  =   051*(C12+C22)-X6AR*(C32+C42+C72+C82>  TRB    171 

Bia6,04)  =   -D51*053*(C12+C22)+XBAR*C;7*(C32+C42+C72+C82)  TRe    172 

8(06,05)  =    (D53+D62)»(C11«(YBAR+C99+C13)+C21*( YBAR-099-022) )  TRB    173 

BIC6.C6)  ■   Cll*( YBAH+D99+D13)**2+C21*( Y8AR-099-C22 )**2+051**2* (C12TRB    174 

*+C22)+XBAR**2*(C32+C42+C72+C82)  TRB    175 

BI06.07)  =    -D51*(  (C62+£14)*C12  +  (C62+C24)»C22)  TRB    176 

S(06,C8)  =    C11*D14*(  YBAR+D99+013)  TRB    177 

EI06.09)  =    C21*C24*(YBAR-099-022)  TR8    178 

B(06,26)  »    C32*XBAR*C34  TRB    179 

£(06,27)  =    C42*XBAR*044  TRB    180 

£(06,28)  =    C72*XEAR*C74  TRB    181 

B(06,29)  "   C82*XBAR*084  TRB    182 

E(07,02)  =    -C12*lQ62  +  D14)-C22*(062+024)  TRB    183 

B(07,03)  =   -C13»(099  +  013)  +  C23*(D99+C22)  TRB    184 

8(07,04)  =    053*(C12*(D62+014)+C22*(C62+024))+YEAR*(C13*(D99+013)-CTR8    185 

*23*IC99+D22>)  TRB    186 

BIC7.C5)  =    D51*(C13*(099+D13)-C23*IC99+C22))  TRB    187 

8(07,06)  =   -051*(C12*(D62+014J+C22*ID62+D24J  )  TRB    188 

3(07,07)  =   C12*(C62+C14)**2+C13*(099+013)**2+C22*( 062+024)**2+C23*TRB    189 

*(099+022)**2  TRB    190 

£(08,01)  =    Cll*014  TRB    191 

B(08,C5I  =    C11*014*(L53+C62)  TRB    192 

E(08,G6)  =    Cll*014*< Y6AR+099+013)  TRB    193 

B(0a,08)  =    C11*C14**2  TRB    194 

B(  09,01)  =   C21*D24  TRB    195 

B(09,05)  =    C21*D24*(C53+C62)  TRB    196 

c(09,06)  =    C21*D24*(YBAR-Q99-D22)  TRB    197 

B109.C9)  =    C21*024**2  TRB    193 

£(10,10)  =    C33  TRB    199 

8(11,11)  =   C43                                                  '  TRB    200 

£112,12)  =   C73  TRB    201 

6(13,13)  =   C83  TRB    202 

6(14,14)  =   C31  TR8    203 

EU4.18)  =    C31*C34  TRB   204 

6(15,15)  =    C41  TRB    205 

6(15,19)  =    C41»044  TRB    206 

£(16,16)  =    C71  TRB    207 

61 16.20)  =  C71*074  TR8  208 
B(17,17)  =  C81  TR8  209 
B(17,21)  =  C81*084  TRB  210 
B( 18,14)  =  C31*C34  TR8  211 
Ella, 18)  =  C31»034**2  TR8  U2 
B( 19,15)  «  C41*C44  TR3  213 
6(19,19)  =  C41*044**2  TRB  214 
8(20,16)  =  C71»C74  TRB  215 
8(20,20)  =  C71*074**2  TRB  216 
6(21,17)  =    C81*D84  TR8    217 

6121.21)  =  C81»084»*2  TRB  218 
B( 26,02)  =   -C32*034  TRB    219 
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E126.04)    =   032*034*057  TRB    220 

8(26, Col    -    C32*D34*XEAR  TRB    221 

E<26,26)    =    032*034**2  TRB    222 

8127,02)    =    -C42*044  TRB    223 

B(27,04)    =    042*044*057  TRB    224 

6127,061    ■    C42*C44»XfiAR  TRB    225 

8(27,27)    =    C42*C44**2  TRB    226 

6(28, C2)    ■    -072*074  TRB    227 

£128,04)    =    C72*C74»057  TRB    228 

8(28,06)    =   C72*074*XEAR  TRB    229 

8(28,28)    =    072*074**2  TRB   230 

B(29,02)    =   -082*084  TRB    231 

8(29,04)    =    082*084*D57  TRB    232 

8(29,06)    =    C82*C84*XEAR  TRB    233 

E(29,29)    =    062*084**2  TRB    234 

C 

0  THE    NEXT    201    CARDS    EVALUATE    ThE    NONZERO    ELEMENTS   OF    THE    MATRIX   C. 

C 

0(01,01)  =  K11+K21+KA1+KA3  TR8  235 

0(01,05)  =  (D53+C62)*(K11+K21)+057*(KA1+KA3)  TRB  236 
C(01,06)  ■  Kll*(YaAR+099+013)+K21*(YBAR-D99-022)+KAl*(056+L53)-KA3TRB  237 

**(D55+(_54)-KA5+KA7  TRB  238 

C101.C8)  =  K11*C14  TRB  239 

CI01.C9)  =  K21*D24  TRB  240 

C(01,14)  =  -KA1  TRB  241 

0(01,15)  =  -KA3        .  TRB  242 

0(01,22)  =  KA1*C33+KA5  TRB  243 

0101,23)  =  -KA3*D42-KA7  TRB  244 

0(02,02)  =  K12+K22+K32+K42+K72+K62  TRB  245 

0(02,04)  =  -D53*(K12+K22)-057*(K32+K42+K72+Ke2)  TRB  246 

0(02, Co)  ■  051*(K12+K22)-XBAR*(K32+K42+K72+K82)  TRB  247 

0(02,07)  =  -K12*(062+014)-K22*(062+024)  TRB  248 

0(02,26)  =  -K32*034  TRB  249 

0(02,27)  =  -M2*D44  TRe  250 

0(02,28)  =  -K72*074  TRB  251 

0(02,29)  =  -K82*084  TRB  252 

C(03,C3)  =  K12+K23+KA1+KA3  TRB  253 

0(03,04)  =  -YdAR*(K13+K23)-(D56+L53)*KAl+(055+L54)*KA3+KA5-KA7  TRB  254 

0(03, C;)  =  -D51*(K13+K23)+XBAR*(KA1+KA3)  TRB  255 

0(03,07)  =  -KU*1099+D13)+K23*(D99+D22)  TRB  256 

0(03,10)  =  -KA1  TRB  257 

0(03,11)  =  -KA3  TRB  258 

0103,26)  =  -KA1*D33-KA5  TRB  259 

0(03,27)  =  KA3*D42+KA7  TRB  260 

0(04,02)  »  -D53*(i<12+K22)-D57*U32  +  K42  +  K72+i<82)  TR8  261 

0(04, U3)  =  -Y8AR*(K13+K23)-(056+L53)*KAl  +  (055+l.54)*KA3+KA5-KA7  TRB  262 
CI  04,04)  =  D53**2*(K12+K22)+YBAR**2*(K13+K23)+057**2*(K32+K42+K72+TR8  263 
*K82)+(056+l53)**2*KAl+(Q55+L54)**2*KA3-2*KA5*{D56+L53) -2*KA7*( 055+TRB  264 

*L54)+KA9+KA11  TRB  265 
C(04,C5)  =  051*Y8AR*(K13+K23)-X6AR*((D56+l53)*KAl-<055+L54)*KA3-KATRB  266 

*5+KA7)  TRB  267 

0(04, C6)  =  -051*D53* (K12+K22 )+XE Afi*C57* (K32+K42+K72+K82 )  TRB  268 
0(04,07)  ■  D53*(  1062  +  014)*K12  +  (062+O24.i-'  Aili  *»aaka*HaS*9-+C13)  *K  13-(  TRB  269 

*099+022)*K23)  TR8  270 

C(04,10)  =  KA1*(D56+L53)-KA5  TR8  271 

C(04,ll)  =  -KA3*(Q55+L54)+KA7  TRB  272 

0(04,26)  ■  034*057*K32+033*(056+L53)*KA1-KA5*(033-C56-L53)-KA9  TRB  273 
C104.27)  =  U44*057*K42+042*(055+L54)*KA3+KA7*(-D42+055+L54)-KA11   TRB  274 

0(04,28)  =  074*057**72  TRB  275 

0104,29)  ■  D84*057*K82  TRB  276 
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CI05,01)    -    (C53  +  C62)»(KU+K21)+D57»(KA1+KA3)  TRB    277 

C(C5,C3)  *  -i)51*(K13+K23)  +  XBAR*(KAl  +  KA3)  TRB  278 
C(05i04J    ■    D51*YBAR*(K13+K23)-XbAR*l ID 56+L53)*KAl- ( D55+L54) *KA3-KATRB   279 

*5+KA7)  TRB  280 
C(05,C5)    ■    <053+D62)**2*(Kll+K21)+L>51**2*  (  K13+K23)  +  (XBAR**2  +  C5  7**2TRB    281 

*)*<KAl+KA3)«KA13+KA15  TRB  282 
CIG5.C6)    =    <D53+D62)*< (YBAR+099+G13 ) *K 11+ ( YB AR-099-C22 ) *K21 ) +D57*( TRB    283 

*(D56+l53)*KA1-(D55+154)*KA2-KA5+KA7)  TRB    284 

C105.07)    =    051*( <099+J13)*K13-(099+C22)*K23)  TRB    285 

C105.C8)    =    014*(D53+C62J*K11  TRB    286 

C105.09)    =    D24*(C52+C6i)*K21  TRB    287 

C(C5,10)    =   -XEAR*KA1  TRB    288 

C(05ill)    ■   -XBAR*KA3  TRB    289 

CI05.14)    ■    -D57»KA1  TRB    290 

C( 05,15)    =   -057*KA3  TRB    291 

C(05,18)    =   -KA13  TRB    292 

C105.19)    =   -KA15  TRB    293 

C(05,22)    =    D57*(KA1*Q33+KA5)  TRB    294 

C105.23)    =    Q57*(-KA7-KA3*D42)  TRB    295 

C(C5>26)    =    -X6AR*D33*KA1-XBAR*KA5  TRB    296 

CI05.27)  =  XBAR*D42*KA3+XBAR*KA7  TRB  297 
CI06.01)    =    Kll*<YBAR+D99+D13)+i<21*<YBAR-0SS-C22)+KAl*<D56+L53>-KA3TRB    298 

**(D55+L54)-KA5+KA7  TRB    299 

C(0e,02)    ■    Q51*IK12+K22)-XBAR*(K32  +  M2  +  K72+K82)  TRB   300 

C(06,04)  =  -051*D53*<K12+K22)+a57*XEAR*(K22+K42+K72+K82)  TRB  301 
C106.C5)    =    (D53+C62)»< ( YBAR+D99+01 3) *K11+ ( YB AR-D99-C22 )*K21 ) +D57*( TRB   302 

*<056+L53)*KAl-(D55+L54)*KA3-KA5+KA7l  TRB  303 
CIC6.C6)     ■    (  YBAR  +  0<,9+013)**2*Kll  +  (YeAR-C99-022)**2»K21+D51**2*(K12TRB    304 

*+K22)+XBAR»*i*(K32+K42+K72+K.a2)  +  lD56+L53)»*2*KAl  +  l055+L54)**2*KA3-TRB    305 

*2*KA5*(C56+L53)-2*KA7*(055+L54)+KA9+KA11  TRB    306 

C106.C7)    =    -D51*(l)62  +  014)*K12-051*<C62  +  D24)*K22  TRB    307 

C106,C8)    ■    C14*IYBAR+D99+013)*KU  TRB    308 

CIC6.C9)     ■    Q24*lYBAR-099-D22)*K21  TRB    309 

CIG6.14)    =    -KA1*(056+L53)  +  KA5  TRB    310 

C106.15)    =    KA3»<G55+1.54)-KA7  TRB    311 

C(06,22)    =    l)33*(D56+L53)  *K A1+KA5*< -C33 +C55+L53 J-KA9  TRB    312 

C106.23)    =    C42»(D55+L54)*KA3+KA7*(-C42+D5c+L54)-KA11  TRB    313 

CIC6.26)     =    X6AR*034*K32  TR8    314 

C(06,27)     ■    XBAR*D44*K42  TRe    315 

C(06,28)    =    XBAR*C74*K72  TRB    316 

C( 06,291    =    X3AR*D34*K82  TR8    317 

CI07.02)    =   -(062+014)*K12-lD62+D24)*K22  TRB    313 

C107.03)  =  -(C99+D13)*K13+(099+022)»K23  TRB  319 
C(07,04)    =    D52»( ( J62+014)*K12+(D62+D24)»K22)+YBAR* ( (099+013 ) *K 13-ITR8    320 

*C99+D22)*K23)  TRB    321 

C(07,05)    ■=    D51*{  (099+D13  )*K13- ( D99+C22 )*K23 )  TRB    322 

C(07,06)  =  -D51*< ( 062+014I*K12+(D62+024)*K22)  TRB  323 
C107.07)    =    (062+C14)»*2*K12+lD99+J13)**2*K13+(D62+024)**2*K22+ ID99TRB   324 

++D22)**2*K23  TRB    325 

C(08,01)    =    Kll*014  TRB    326 

C(08,05)    «    Kll*01"t*(C53+C62>  TRB    327 

0(08,06)    ■    K11*D14*( YBAR+D99+013)  TRB    328 

C(03,06)    =    <11*D14**2  TRB    329 

K21*024  TRB    330 

K21*024*(D53+062)  TRB   331 

K2l*C24*(YBAR-099-022)  TRB    332 

K21*024**2  TRB    333 

-KA1  TRB   334 

KA1»(C56+L53)-KA5  TRB    335 

-KA1*XBAR  TRB    336 
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C( 10,10)     =    K33+KA1+KA2  TRB    337 

CI  10,  12)    =    -KA2  TRB    338 

C(10, 26)    =    KA1*C33-KA2»(C32-H.37I*KA5+KA6  TRB    339 

CUC, 28)    ■   -KAi*Q73-KA6  TRS    ^° 

CI11.03)    =    -KA3  TRB    341 

C(H, 04)    =   -KA3*(055*L54)+KA7 

CI  11,05)    =    -KA3*XBAR 

C(ll.ll)    =    K43+KA3+KA4 

CU1.13)     •    -KA4 

CI  11,27)    =    -KA3*D42+KA4*(D43+L84)-KA7-KA3' 

C111.29)    =    KA4*C82+KA8 

0(12,10)     =    -KA2 

C(12,12)    =    K73+KA2 

0(12,26)    =    KA2*(C32  +  L37)-KA6 

C( 12,28)    =    KA2*D73+KA6 

CI  13, 11)    =    -K.A4 

C113.13)     =    K83+KA4      ' 

CI  13,271    =    -K.A**043+Le4)+KA8 

CH3.29)    =    -KA4»D82-KA8 

C( 14,01)     =    -KA1 

0(14,05)    =    -KA1*D57 

C(14,06)    ■    -KA1*(056+L53)+KA5 

C( 14,14)    =    K31  +  KA1+KA2 

CI14.16)    =    -KA2 

C(14,18)     =    K31*034 

C(  14,22)    ■   -K.A1*033+KA2*(D32K.37)-KA5-KA6 

C114.24)    =    +KA2»C73+KA6 

C(15,01)     =    -KA3 

C(15,05)    =    -KA3*D57 

C(15,06)    =    KA3»(C55  +  L54)-KA7 

C(  15,15)     =    K4H-KA3+KA4 

C115.17)    =    -KA4 

C( 15,19)     =    K41*044 

0(15,23)    =    K.A3*C42-KA4*(  D43+L84) +KA7+KA8 

0(15,25)    =    -KA4»C82-KA8 

C(  16,14)     =    -KA2 

C(16,lt)    =    K71+KA2 

C(16,20)     =    K71*D74 

C(lt, 22)    *  -KA2*(J32+L37)+KA6 

0116,24)    =    -KA2*C73-KA6 

C(17,15)     =    -K.A4 

0(17,17)    =    K81+KA4 

0117,21)    =    K81*C84 

0(17,23)    =    KA<»*  (043-1- 184) -KA8 

C( 17,25)    =    KA4*C82*KA8 

0(18,05)    =    -KA13 

C(18, 14)    =    K31»034 

0(18,18)    =    K.31*C34**2+KA13+KA14 

C( 18,20)    =   -KA14 

C(19,05)    =    -KA15 

C(1S,15)     =    K41*C44 

C(  19,19)     =    K41*044»*2+KA15+KA16 

C(19,21)  =  -KA16 

C(20,l6)  '■  K7l*074 

C(20, 18)  =  -KA14 

0120,20)    =    K71»074*»2+KA14 

CI21.17)    =    Kdl*084 

0(21,19)    =    -KA16 

C(21,21)     •    K61*C84**2+KA16 

0(22,01)  =  KA1*033+KA5 
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C(22,05)  =    KA1*C33*057+057*KA5  TRB    397 

0(22,06)  ■    (D56+L53)*(KA1*D33<-KA5)-KA5*D33-KA9  TRB    398 

0(22,14)  =    -KAl»C33+KA2*(032+i.37)-KA5-KA6  TRB    399 

0(22,16)  =   -KA2*(032+L37)+KA6  TRB    400 

0(22,22)  =    KA1»033**2+KA2*(D32+L37)**2  +  2*KA5*D33-2*KA6*(C32+L37)+KTRB   401 

♦A9+KA10  TRB    402 

C(22t24)  »    KA2*D73*(032+L37)+KA6*(-073-»032+L37)-KA10  TRB    403 

C(23,01)  =   -KA3*042-KA7  TRB   404 

C(23,05)  =   -KA3*C42*C5  7-C57*KA7  TRB    405 

0(23,06)  =    (D55-H.54)*!  KA3*D42+KA7) -KA7*D42-KA1 1  TRB    406 

0(23,15)  =    KA3»D42-KA4*(C43+L&4)+KA7-tKA8  TRB    407 

0(23,17)  =    KA4«(D43+L84)-KA8  TRB    408 

0(23,23)  =    KA3*D42**2+KA4*<D43*L84)»*2+2*KA7*042-2*KA8*(D43+L84)+KTRB    409 

♦A11+KA12  TRB    410 

C123.25J  =    KA4*D82*(D43+L84)-KA8*(CE2-C43-L84)-KA12  TRB    411 

0(24,14)  =    KA2*C73+KA6  TRB   412 

0124,16)  =    -KA2*D73-KA6  TRB    413 

C(24,22)  ■    (D32+L37)*(KA2*073+KA6)-KA6*073-KA10  TRB    414 

C(24,24)  =    KA2*C73**2+2*KA6*D73+KA10  TRB    415 

0(25,15)  »  -KA4*G32-KA8  TRB    416 

0(25,17)  =    KA4*082+KA8  TRB    417 

0(25,23)  =     (D43+L34J*(KA4*C82+KAa)-KA3*C82-KA12  TRB    418 

C( 25,25)  =    KA4*D82**2+2*KA8*082+KA12  TRB    419 

C(26t02)  ■    -K32*C34  TRB    420 

C(2a,03)  =   -KA1*D33-KA5  TRB    421 

C(26,04)  ■    K32*057*034  +  (056+(.;3)*(KA1*033  +  KA5)-KA5*C33-KA9  TRB    422 

C(26,C5I  =    -KA1»C33*XBAR-XEAR*KA5  TRB    423 

C(2o,06)  =    K32*Q34*XBAR  TRB    424 

C(26,10l  =    KA1*C33-KA2*(032+L37)+KA5*KA6  TRB    425 

0126,12)  =    KA2»(C32+L3  7)-KA6  TRB    426 

0(26,26)  =    K32*034**2+KAl*D33**2+KA2*<032-i-L37)**2+2*KA5*D33-2»KA6»TRB    427 

*  (D32-fL37)+KA9+KA10  TRB    428 

0(26,26)  ■    KA2*(D32+L37)*073-KA6*< C73-C32-L37 ) -KA10  TRB    429 

C(27,02)  =    -K42*D44  TRB    430 

C127.03)  =    KA3*C42+KA7  TRB    431 

C127.04)  =    K42*(344*D57+(055+1.54J*(  KA3*D42+KA7) -KA7*Q42-KA11  TRB    432 

C(27,05)  =    KA3*C42*XEAR+XBAR*KA7  TRB    433 

C(27,Col  »    K42*D44*XEAfl  TRB    434 

C(27,ll)  ■   -*A3*D42+KA4*(043<-L84)-KA7-KA8  TRB    435 

0(27,13)  -    -KA4*(J43+L64)  +  KA8  TRB    436 

C(  27,27)  =    K42*D44**2+KA3*042**2+KA4*( 043 +L84) **2+2*KA7»042-2* KA8*TRB    437 

*(043+L84)+KAll+KA12  TRB    438 

0(27,29)  ■    KA4*((J43  +  LS4)»082-KA8*(Ce2-043-Ld4)-KA12  TRB    439 

C(28,02)  =    -K72*074  TRB    440 

0(28,04)  =    K72*C74*D57  TRB    441 

C( 28,06)  =    K72*074*XEAR  TR8    442 

C128.10)  =    -K.A2*D73-KA6  TRB    443 

C(28,12)  =    KA2*073+KA6  TRB    444 

C(28,26)  =    KA2*073*(G32+L37)+KA6*(032+l.37-073)-KA10  TRB    445 

C(28,28)  ■    K.72*C74**2+KA2*D73**2+2*KA6*073+KA10  TRB    446 

0(29,02)  =    -K62»Ca4  TRB    447 

C(29,04)  =    K82*084*D57  TRB    448 

0(29,06)  =    :;82*C84*XEAR  TRB    449 

0(29,11)  =    KA4*082+KA8  TRB    450 

0(29,13)  =    -KA4*C32-KA8  TRB    451 

0(29,27)  =    KA4*ca2*(C43«-U>4)  +KA8*< D43+L34-082J-KA12  TRB    452 

0(29,29)  =    K82*084**2+KA4*082**2+2*KA8*082+KA12  TRB    453 

C  THE    NEXT  113    CARDS    EVALUATE    Tt-.E    NONZERO    ELEMENTS    OF   THE   MATRIX    D. 

C 
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0(01,01)  =  -1.  TRB  454 

0(01,04)  =  -Kll  TRB  455 

C(01,07)  =  -K21  TRB  456 

0(01,22)  "  -Cll  TRB  457 

0(01,25)  =  -C21  TRB  458 

0(02,02)  =  -1.  TRB  459 

0(02,05)  =  -K12  TRB  460 

C(O2,0B)  =  -K22  TRB  461 

0(02,11)  =  -K32  TRB  462 

C(02,14)  =  -K42  TRB  463 

C102.17)  =  -K72  •  TRB  464 

0(02,20)  =  -K62  TRB  465 

0(02,23)  =  -C12  TRB  466 

0(02,26)  ■  -C22  TRB  467 

0(02,29)  ■  -C32  TRB  468 

0(02,32)  =  -C42  •  TRB  469 

0(02,35)  =   -C72  TRB  470 

0(02,38)  =  -C82  TRB  471 

0(03,03)  =  -1.  TRB  472 

01 03, C6)  ■  -K13  TRB  473 

0(03, 09)  =  -K23  TRB  474 

0(03,24)  =  -C13  TRB  475 

0(03,27)  =  -C23  TRB  476 

0(04,02)  =  CB3  TRa  477 

0(04,03)  =  YEAR  TRB  478 

0(04, C5)  =  K12*D53  TRB  479 

0(04,06)  =  K13*YEAR  TRB  480 

D(04,Cd)  =  K22*053  TRB  481 

0(04, C9I  =    K23*YBAR  TRB  482 

0(04,11)  =  K32»C57  TRB  483 

0(04,14)  =  K42»U57  TRB  484 

0(04,17)  =  K72*C57  TRB  485 

0(04,20)  ■  K62*057  TRB  486 

0(04,23)  =  C12*053  TRB  487 

C104.24)  =  C13*YBAR  TRB  488 

0(04,26)  =  C22*053  TR8  489 

0(04,27)  =  C23*YBAR  TRB  490 

0(04,29)  =  C32*D57  TR3  491 

0(04,32)  ■  042*057  TRB  492 

0(04,35)  =  C72*C57  TRB  493 

0104,38)  =  082*057  TRB  494 

0(05,01)  ■  -0B3  TRB  495 

0(05  ,03)  ■  -061  TRB  496 

0(05,04)  =  -K11»(D53+D62)  TRB  497 

0(05,06)  =  K13*051  TRB  498 

0(05, C7)  =  -K21* (053+062)  TRB  499 

0(05, C9)  =  K23*D51  TRB  500 

0(05,22)  -  -l_ll*(053*D62)  TRB  501 

0(05,24)  =  C13*05l  TRB  502 

0(05,25)  =  -C21*(053+D62)  TRB  503 

0(05,27)  =  C23*D51  TRB  504 

0(06, CI)  =  -Y3AR  TRB  505 

0106,02)  =  081 

0(06, C4)  =  -K11*(YBAR+C99+D13) 

0(06,05)  =  -K12*051 

0(06,07)  ■  -K21*(YEAR-C99-D22J 

D(06,Ca)  •  -K22*D51 

0(06,11)  =  K32*XEAR 

0(06,14)  =  K42*XEAR 

0(06,17)  «  K72*XBAR 


TRB 
TRB 

506 
507 

TRB 

508 

TRB 

509 

TRB 

510 

TRB 

511 

TRB 

512 

TRB 

513 

97 


0(06,20)     =    K82*XEAR 

0(06,22)    =    -C11*(YBAR+D99+D13) 

0106,23)    =    -C12*D51 

0(06,25)    «   -C21*(YBAR-C99-D22J 

0(06,26)    ■    -022*051 

0(06,29)    =    C32*XEAR 

0(06,32)    =    C42*XBAR 

0(06,35)    =    C72»XEAR 

D(06,38)     ■    C82*XEAR 

0(07, C5)    =    K12*(D62+014) 

C107.C6)    =    K13*(D99+C13) 

D107.C6)    ■    K22*(D62<-024) 

0(07,09)    =   -K23*(D99*D22) 

0(07,23)    =    CI  2* (062+ 014) 

0(07,24)    =    C13*(D99+013) 

0(07,26)    =    C22*(C62*024I 

0107,27)    =  -C2j*(D99+022J 

0(08,04)    =    -K11*D14 

0(08,22)    ■    -Cll*014 

0(09,07)    =   -K21*024 

0(09,25)    =    -021*024 

0(10,12)     =    -K33 

0(10,30)    =   -C33 

0(11,15)    =    -K43 

D( 11,33)    =   -C43 

0(12,18)    =   -K73 

0(12,36)    *    -C73 

D(  13,21)    =   -K83 

01 13,39)    =    -C83 

0(14,10)     =    -K31 

01  14,23)    =   -C31 

D(15,13)    =    -K41 

U(15,31)     =    -C41 

0(16,16)    ■    -K71 

0(16,34)     =    -C71 

0(  17,  19)    ■    -K81 

0(17,37)    =    -C81 

0(16,10)    ■    -K31*D34 

0118,26)    =    -C31*D34 

0(19,13)    =    -K41*C44 

D( 19,31)    =    -041*044 

0(20,16)    =    -K71*D74 

DI20.34)    ■    -071*074 

0(21, 19)    =    -KE1*084 

0121,37)    =    -081*084 

0(26,11)    =    K32*034 

0(26,29)    ■    032*034 

0(27,14)    ■    K42*C44 

DI27.-2)    =    042*044 

D(28,171    =    K72*074 

0(28,35)    *    072*074 

D(29.2C)    »    K82*0£4 

0(29,38)    =    082*084 

c 

c 
c 

THE    NEXT    10   CARDS   CAST    THE    A, 

00    60    1*1,29 

DO    40    J=l,29                   i 

Gil, J)    =    -A1I.JI 

Gll.J+29)    =    CH.J) 

TRB  514 
TRB  515 
TRB  516 
TRB  517 
TRB  518 
TRB  519 
TRB  520 
TRB  521 
TRB  522 
TRB  523 
TRB  524 
TRB  525 
TRB  526 
TRB  527 
TRB  528 
TRB  529 
TRB  530 
TRB  531 
TRB  532 
TRB  533 
TRB  534 
TRB  535 
TRB  536 
TRB  537 
TRB  538 


TRB    539 

TRB    540 

TRB    541 

TRB    542 

TRB    543 

TRB    544 

TRB    545 

TRB    546 

TRB    547 

TRB    548 

TRB    549 

TRB    550 

TRB    551 

TRB    552 

TRB    553 

TRB    554 

TRB    555 

TRB    556 

TRB    557 

TRB    558 

TRB    559 

TRB    560 

TRB    561 

TRB    562 

TRB    563 

TRB    564 

TRB    565 

TRB    566 

HE    MATRIX    G. 

TRB    567 

TRB    568 

TRB    569 

TRB    570 
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GlI.J+58)  =  BII.JI  TRB  571 

40  CONTINUE  TR8  572 

CO  50  J=l,39  TRB  573 

G(I,J+87)  =  DU,J)  TRB  574 

50  CONTINUE  TRB  575 

60  CONTINUE  TRB  576 

CALL  GAUSS  (G.N.Nl)  TRB  577 

WRITE  (10)  G  TRB  578 

C 

C  THE    NEXT    13    CARDS    FRCOUCE    A    FORMATTEC    COPY    OF   THE   NCN^ERO    ELEMENTS 

C      OF  TFE  MATRIX  G. 

C 

DO  80  1=1,29  TRB  579 

WRITE  (6,1501  I  TRB  580 

CO  70  J=l,126  TRB  581 

IF  U.ea.il  ID  =  1  TRB  582 

IF  (ABS(G(  I, J)).LE. 0.0001)  GO  TO  70                                TRB  583 

IF  (ID.EC.l)  MUTE  (6,90)  G(I,J),NAME(J)  TRB  584 

IF  (ID.EQ.2)  WRITE  (6,100)  G(  I  ,  J  )  ,  NAME  (  J)                          TRB  585 

IF  (ID. E'J. 3)  WRITE  (6,110)  G(  I ,  J  ) ,  NAME  U  )                          TRB  586 

IF  (ID.E0.4)  WRITE  (6,120)  G ( I , J) , \AKE ( JJ                          TR8  587 

IF  (IC.EQ.4)  ID  =  1  TRB  588 


ID  =  ID+l 


TRB  565 


70  CONTINUE  TRB  590 

80  CGNTINUE  TRB  591 

STOP  TRB  592 

90  FORMAT  ('  ' , 20X.F3.1 , IX, A4, < =' )  TRB  593 

100  FORMAT  I»+',T3l,'+('  ,F13.4,'  )',A4)  TRB  594 

110  FORMAT  (•+' ,T51,'+C  .F13.4,' )' ,A4)  TRB  595 

120  FORMAT  ( ' + • ,T71 , '+( ■ ,F 13 .4 , ■ ) ■ ,A4, /)  TRB  596 

130  FCRMAT  (20A4)  TRB  597 

140  FORMAT  (8F10.2)  TRB  598 

150  FORMAT  ( '-' , 47X, • ECUAT ION ' , I  5)  TRB  599 

160  FCRMAT  (2E10.4)  TRB  600 

ENO  TRB  601 


99 


SUBROUTINE   GAUSS    (A.N.M)  GJR      10 

C  THIS    PROGRAM    FEPFQRMS    A    GAUSS-JGROAN   REDUCTION    ON    AN   N    X  Nl    MATRIX   A. 


C 


DIMENSION    A(N,N1) 

CO    SO    J=1,N 

£IG    •    ABSIAIJ.JJ) 

KSAVE    =    J 

IF    (J.EC.N)    GC    TO    30 

KK   '    J  +  l 

00    10    K=KK,N 

IF    (AeSU(K,J)).LE.BIGl    GO   TO    10 

BIG    ■    ABS(A(K,J)J 

KSAVE    =    K 
10    CONTINUE 

00    20    f.=  J,Nl 

DUMMY    =    A(KSAVE,M) 

AIKSAVE.M)    =    AiJ.M) 
20    A(J.M)     =    OUMMY 
30    OIV    =    A(J,J) 

S    =    1.0/DIV 

00    40    K=J,N1 
40    A(J,KJ    ■    A(J,KI*S 

00    70    1=1, N 

IF    (I-JJ    5O.7C50 
50    AIJ   =    -A (I, J) 

DO    60    K=J,N1 
oO    Ad.KJ    =    A( I,K)+AIJ*A( J,K) 
70    CONTINUE 
30    CONTINUE 

RETURN 

END 


GJR 

20 

GJR 

30 

GJR 

40 

GJR 

50 

GJR 

60 

GJR 

70 

GJR 

80 

GJR 

90 

GJR 

100 

GJR 

110 

GJR 
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GJR 
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160 
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270 

GJR 

280 
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THIS    PROGRAM    INTEGRATES    The    EQUATIONS    OF    MGTION    USING    THE    RUNGE-KUTTA- 
GILL    METHOD. 

DIMENSION    G<29,126),    HI58,97),    g(58J,    XSE),    DX(58),    Y(  39)  INT         1 

CGMMON   H,Y  INT         2 

REWIND    10  INT         3 

THE    MATRIX   G    IS    READ    FROM    OISK. 

READ    (10)    G  INT         4 

THE    NEXT    THREE    CARDS    READ    IN    ThE    INITIAL    CONDITIONS,    INCREMENT    SUE, 
AND    NUMBER    OF    STEPS    TO    BE    TAKEN. 

REAC  (5,90)  (XII),  1  =  1,58) 

READ  (5,50)  DELT 
READ  (5,100)  NSTEP 

THE  NEXT  4  CARDS  INITIALIZE  THE  H  MATRIX. 

CO    20    1=1,58 

DO  10  J=l,97 
10  h(  I, J)  ■  0.0 
20    CONTINUE 

ThE   NEXT    5    CAROS    CAST    THE    G    MATRIX    INTO    ThE    H    MATRIX.    H    CONTAINS    THE 
EXPRESSIONS    FOR    THE    DERIVATIVE    OF    EACH    X. 

DO    40    1=1,29 
H(i ,1+29)     «    1.0 
DO    30    J=l,57 
30    H(I+29,J)    =    GII.J+29) 

40    CONTINUE 

THE    NEXT    4    CAROS    INITIALIZE    THE    WORK    VECTOR    Q,    AND    THE     SURFACE 
CONDITION    VECTOR    Y. 

DO    50     1=1,53  '                                                                                                                                  INT  17 

50  am    =   0.0  INT  18 

DO    60    1=1,39  INT  19 

60    Y(I)     =   0.0  INT  20 

NEQ    IS    ThE    NUMBER    OF    FIRST    ORDER    EQUATIONS    TC    BE    INTEGRATED. 

NEQ   ■'   58  INT      21 

THE    TIKE    IS    INITIALLY    0.0    . 

T    =    0.0  INT      22 


INT 

5 

INT 

6 

INT 

7 

INT 

8 

INT 

9 

INT 

10 

INT 

11 

INT 

12 

INT 

13 

INT 

14 

INT 

15 

INT 

16 

c 

c 

c 
c 
c 

c 

C  THE    NEXT    4    CARDS    PRODUCE    A    HARD    COPY    OF    THE    INITIAL    TIME    AND    DISPLACEMENTS 

C  ,    ANO    WRITE    THE    VARIABLE    T    AND    VECTOR    X    ON    DISK. 

C 

WRITE    (6,120)    T 

WRITE     (6,110)     (XII), 1=1, 29) 

WRITE    (10)    T 

WRITE    (10)    X 
C 

C  THE    NEXT    8    CARDS    CALL    THE    RKG    SUBROUTINE.       THE    PROCEDURE    IS     PERFCRMEO 

C  TEN    TIMES    FOR    EVERY    STEP    REQUESTED    BY    NSTEP.       T    ANO    X    ARE    PRINTED 


INT 

23 

INT 

24 

INT 

25 

INT 

26 

0OT  AND  STGREC  CN  DISK. 
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DO   80    J=1.NSTEP 
DO    70    N=lilO 

CALL    RKG    (NEQ.D£LT,T)X,OXlQ) 
70   CONTINUE 

WRITE    (10)    T 
WRITE    (10)    X 
WRITE    (6,120)    T 
80    WRITE    (6,110)    (X(I),  1=1,29) 

STOP 
90    FORMAT    (6F10.4) 
100    FORMAT    (13) 
110   FORMAT    (•     « .12F10.4) 


120    FORMAT 
END 


1 '-'.F13.5I 


INT 

27 

INT 

28 

INT 

29 

INT 

30 

INT 

31 

INT 

32 

INT 

33 

INT 

3* 

INT 

35 

INT 

36 

INT 

37 

INT 

38 

INT 

39 

INT 

40 
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SUBRGLTINE    «KG<NEy,H,X,Y,UY, S)  RKG       10 

C  ThE    INDEPENDENT    VARIABLE    X    IS    INCREMENTED    IN    THIS    PROGRAM 

C  Yd)     AND    DY(I)     ARE    TFt    DEPENDENT    VARIABLE    AND    ITS    DERIVATIVE 

C  ALL    THE    Q(I)    MLST    BE    INITIALLY    SET    TC    ZERC    IN    THE    MAIN    PROGRAM 

C  h  ■    INTERVAL    SIZE 

C  NEQ    =    NUMBER    CF    FIRST    CRCER    EGUATICNS 

C  A    SUBROUTINE    DERIVINEQ ,X , Y, DY)    MLST    BE    PROVIDED 

C 

DIMENSION      A(2) 

DIMENSION    Y(NEQ),DY(NES) ,0(NEQ) 

A(l)=C.292aS3216  813452475 

A(2)=1.7C71C6 7811865 47 52 

H2=.5*h 

CALL    DERIV(NEC,X,Y,DY) 

DO    13    1=1, NEQ 

B=H2*0Y(I)-g(I) 

Yl  I)=Y(I)+B 
13    Qd)  =  Qd)+3.*B-h2*CYd) 

X=X+H2 

CO    20    J= 1,2 

CALL    CERIV(NEC,X,Y,0Y1 

CO    20    1=1, NEQ 

E=A(J)*(h*DY(I)-Q(HJ 

Yd)  =Y(I)+B 
20    Q(I)=C(I)+3.*B-A(J)*h*0YII) 

X=X+H2 

CALL    DERIV(NEC,X,Y,DY) 

DO    26    1=1, NEU. 

E=0. 106666666 666666666* I H*DY ( I )-2.*C ( I ) ) 

Y(  I)  =  Yd)+B 
26    Q( I)=C(I)+3.*B-H2»0Y(I) 

RETURN 

END 
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SUBROUTINE  DERIVCNEQ.T.X.OX)  OER   10 

C      THIS  SUBROUTINE  EVALUATES  ThE  DERIVATIVE  OF  EACH  X  FOR  USE  IN  THE  RKG 

C      SUBROUTINE. 

C 

CIMENSICN  X(NE0),DX(NEO  ,H(58,97),  Y<39)  OER   20 

COMMON  H,Y  OER   30 

C 

C  THE   NEXT    11    CARCS    EVALUATE    THE    SURFACE    CONOUTIGNS   FOR    THE    LEFT    REAR 

C  WHEELS    TRAVERSING    A    SINUSCICAL    BUMP. 

C 

IFIT.GE. 0.6815)    CO    TO    10 

Y(12)=-0.416*SIN(4.61*T) 

Y(18)=Y(12) 

Y(30)=-1.918*COS(4.61*T.) 

Y(36)  =  Y(30) 

GO    TO    20 
10    Y(12)=0.0 

Y(18)=0.0 

Y130)=0.0 

Y(36)  =  0.0 
20    CONTINUE 


THE    NEXT    7   CARDS    EVALUATE    THE    DERIVATIVES. 

00  50  1=1,58 

CX(I)=0.0 

DO  30  J=l,58 
30  DX(II=DX(IJ+H(I ,JJ*X(J) 

DO  40  J=l,39 
40  0X(1)«0X(I)+H(I ,J+58)*Y(J) 
50  CONTINUE 

RETURN 

ENO 
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APPENDIX  B 


DERIVATION  OF  THE  POTENTIAL  ENERGY  IN  THE  REAR  AXLES 
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The  rear  axles  are  treated  as  cantilever  beams  with  lumped  masses  at 
their  ends,  as  in  Figure  Bl 


FIGURE  Bl 


The  lumped  mass  causes  a  force  F  and  moment  M  to  be  transferred  to  the  end  of 
the  beam.   The  deflection,  x,  and  slope,  8,  at  the  end  of  the  beam  are 
obtained  using  elementary  beam  theory  and  superposition. 


x  «  FL3/3EI  +  ML2/2EI 
8  =  FL2/2EI  +  ML/EI 


(Bl) 
(B2) 


Equations  Bl  and  B2  may  be  used  to  find  expressions  for  F  and  M  in  terms 
of  x  and  8. 


F  =  12EI/L3  x  -  6EI/L2  8 

M  -  -  6EI/L2  x  +  4EI/L  8 

The  differential  of  energy  for  the  beam  is  given  by 

dU  "  Fdx  +  Md8 


(B3) 
(B4) 

(B5) 
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dF/d9  is  equal  to  dM/dx,  therefore  B5  is  an  exact  differential,  which 
implies 

dU/dx  =  F  (B6) 

and 

dU/de  =  M  (B7) 

Substitution  of  (B3)  into  (B6)  and  integrating  yields 

U  =  6EI/L3  •  x2  -  6EI/L2  •  x-9  +  f(6)  (B8) 

Substituting  (B8)  into  (B7)  yields 

dU/d8  -  -  6EI/L2-x  +  f'-(9)  (B9) 

which  must  correspond  to  (B4).  To  do  this, 

f'(e)  =  4EI/L-9  (BIO) 

Integrating  (BIO)  gives 

f(8)  =  2EI/L-92  +  C  (Bll) 

where  the  constant  C  is  a  reference  level  which  may  be  set  to  zero  for  this 
problem.  Substituting  (Bll)  into  (B8)  yields  the  final  expression  for  the 
energy  in  the  axle 

U  -  6EI/L3-x2  -  6EI/L2-x6  +  2EI/L-92  (B12) 
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APPENDIX  C 


SAMPLE  EQUATIONS  OF  MOTION  AND  INTEGRATION  RESULTS 
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EQUATION  1 

1.0   AX1   =*(      -29185.66411X1      ♦(  969.56591X3      »(  - 17305. 6992) X5 

+  (           -128.53171X8      +  <  -126.53171X9      +{  -492. 249d) X 10 

+(            -492.24931X11    ♦  <  14478. C7C3 ) X14    +(  14476.07031X15 

-50.96221X18    ♦(  -50.96221X19   *<  - 19255.  <S242    X22 

+  <          19255.82421X23    *■<  -654.69241X26    +<  654.69241X27 

+  <                 -1.26241VX1    +(  -0.12291VX3    ♦  <  -1.5056IVX5 

+  (                -0.7C69IVX8   +(  -0.706S1VX9   +{  O.00731FCB1 


+(  -0.0003)FQB3+(  114.76051Y11    +(  7 

+(  114.76051Y21    +1  7.46741Y23    +( 

+(  0.0615)YD13+(  0.63121Y021+1 


4674)  Y13 
0.63121 YOl 1 
0.06151YD23 


EQUATION  2 

1.0    AX2    .+{  -348.98561X2      ♦(  -35122.52341X4      ♦«       -16947.16801X6 

*  140.6S821X7      +(  -12.48411X8      +  <  12.48411X9 
t  "J?!?*«»f J5J2  *!  «27. 26561X11   +(           2163.65011X14 
+  -2l6j. 85011X15    +(  -2689.59591X22+1         -2669.5  9591X23 
+  -5513.56251X26    H  -5515.56251X27   +(              111.66131X28 

♦  111. 68131X29    ♦(  -C.89421VX2    ♦  <  0.8574    VX4 
+  0.02071VX6    ♦  <  C.32671VX7   +  <                -0.0687    VX8 
+  0.06671VX9   ♦  <  C.3CC3)VX26+(                   0.30C3    VX27 
+(  0.3003)VX28+(  C.30C31VX29+1 
+  <  11.14661Y11    ♦  (  62.8U71Y12    +( 
+  (  62.8117JY22   +(  55.S4C61Y32   ♦  ( 
+  (  55.8406JY72    ♦  <  55.84061Y82   +  ( 

0*JSSSiI°"+'  -0.0613)Y021+(                    0.146e)YD22 


0.  00461FD82 
-11.1466)  Y21 
55.  8406)  Y42 


+  <  0.1502)YD32+<  C.  1502  )  YD42K 

+<  0.1502JYD82 


.    .   ...  EQUATION  3 

1.0  AX3    =+(  -529.21831X1      +(      -31801.37501X3      +( 

+1  -5.92101XO      +(  -5.92101X9      <■  ( 


0.1502) YD72 


72931.56251X5 
15753. 1133) X10 


+(  15753.11381X11   ♦<  259.3223JX14  ♦(  259.32231X15 

♦j  120.71281X18    +(  120.71281X19    ♦<  -344.  6S89)  X22 

t  n"o==?lX23    *'  20951.63281X26    ♦(  -20951.63281X27 

+  -0.05821VX1    +  <  -2.429DVX3    ♦{  8.82451VX5 

*i  -0.03261VX8    *J  -0.03261VX9    ♦< 


0.0C0DFCB1 


♦J  0.0C79)FDB3+<  S.iUuillSl  147.' 5721 )  YU 

+  (  5.28661Y21    +1  147.572DY23    +( 

+  1  1.2145JYD13-M  0.02911Y021-M 


0.  0291) YD1 1 
1.2145)  Y023 


EQUATION  4 

t                 li,VlihrZ\*Z      *'    -1""7-J125)X4  *(  472.361DX6 

+  <               -77.8404)X7      +  1                    0.244o)X8  +(  -0.24461X9 

+         "18830.12891X10    +(          18830.1289  1X11  +(  -42.40431X14 

+                  42.40431X15   ♦ (                52.70721X22  +(  52.70721X23 

♦  (       -23430.12501X26    M       -23430.1250  1X27  +<  -llAilVnZi 
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EBUATION  9 

1.0    AX9    =+(       -15719.30081X1      +  (      -24364. 9023 ) X5      +( 
+  (       -17605.60941X9       +  (  -86.45621VX1    +  < 

+  (  212.68211VX6    +  (  -96.83091VX5    +( 

+  (  86.45621YC21 


38669. 468a)X6 

-134.0070)VX5 
15719. 3008)Y21 


EQUATION         10 
1.0    AX10=+(       176081.12501X3      +  <    -674390.8125  1X4      +( 
+  <    -595268.43751X10    * (       4173 77. £7 50 1X12   + ( 
+(      450767.61251X28    +<  -21.8888 1VX10+ ( 

+(  21.8883IYC33 


410269.06251X5 
216580.25001X26 
1309.35551Y33 


EQUATION  11 

1.0  AXll=+l  176081.12501X3  +  (  o74390. 8125 1X4  +  ( 
+  (  -595268.43751X11  +  (  417377.67501X13  +  < 
+  (    -450767.81251X29    +  <  -21. 3888 1VX11+ ( 

+(  21.88881YC43 


410269.06251X5 
-216530.25001X27 
1809.35551Y43 


EQUATION         12 
1.0    AX12=+(       417377.87501X10    +  (    -419187.2500  1X12    +  < 
+  <    -450767.O1251X28    + (  -21. 38 88 1 VX 12+ ( 

+(  21.68831YC73 


-450768.25001X26 
1809.35551Y73 


EQUATIGN  13 

1.0    AX13=+<       417377.87501X11    +(    -419187.25001X13    +( 
+(      450767.81251X29    +(  -21. 86 88 1VX13+ ( 

+  (  21.88881YD83 


450768.25001X27 
1809.35551Y83 


EQUATION  14 

1.0    AX14=+(       176081.12501X1       +(       117974 .25C0 ) X5       +( 

+(    -595047.75001X14    +(       4173  77.37501X16    +( 

+(    -216530. 25001X22    +(    -450767.3125  1X24   +( 

+(  -23.65401VX18+ (  1588. 7024 1 Y31    + ( 


674390.81251X6 
-3177.40481X18 
-11.82701VX14 
11.82701YD31 


EQUATION         15 

.0    AX15=+(       176081.12501X1       +(       117974.2500  1X5  +( 

+1    -595047.75001X15    +(       417377.67501X17  +( 

+(       216580.25001X23   +(      450767.6125  1X25  +( 

+(  -23.65401VX19+I  1588 . 7024 1Y4 1  +{ 


-6 74390. 8125JX6 
-3177.40481X19 
-11.82701VX15 
11.82701YD41 


EQUATION  16 

.0    AX16=+(        417377.67501X14    +(    -418566.56  251X16  +( 

+(       450768.25001X22    +(       450767.81251X24  +( 

+(  -23.6540IVX2C+1  1588. 7C 24 1 Y7 1  +1 


-3177.40461X20 
-11.  82701VX16 
11.82701Y071 


EQUATION  17 

.0  AX17=+(  417377.87501X15  +(  -418566.56251X17  +( 
+(  -450768.25001X23  +1  -450767.81251X25  +1 
+(  -23.6540JVX21+I  1586. 7C24) Y31    +( 


-3177.4048)  X21 
-11.8270)  VX17 
11.  8270)  YD81 


109 


+  (  -24.38791X29    +  1  0.2S621VX2    +  (  -0.34141VX4 

+  (  0.2900)VX6    +  <  -0.1819IVX7    +  (  0.00131VX8 

+(  -0.00131VX9    +(  -0.06691VX26+1  -0.06691VX27 

+(  -0.06691VX28+I  -0.06691VX29+I  -0.00101FC82 

+(  -0.21841Y11    +(  -34.75021Y12    +(  0.21841Y21 

+(  -34.75021Y22    + (  -12.44391Y32    +  <  -12.44391Y42 

+(  -12.44391Y72   +  <  -12.44391Y82    +(  -0. 0012 ) YD1 1 

+  (  -0.08121YC12+I  0.0012)Y021+(  -0.08121YD22 

+(  -0.0335)YC32+(  -0. 0335 1 YD42+ (  -0.0335)YD72 

+(  -0.0335IYC82 


EQUATION  5 

1.0    AX5    =+<  -2251.25271X1      +  <  -9765.6523  1X3      +  (  -26231.39451X5 

+  <  -25.18761X8      +  (  -25.18761X9      +  <  4958.03521X10 

+  (  4958.03521X11+1  1103.13701X14+1  1103.13701X15 

+(  513.50241X18    +(  513.50241X19    +(  -1467.17241X22 

+(  1467.17241X23+1  6594. 16 75 1 X2e    + (  -6594.18751X27 

+(  -0.24741VX1    +(  1.23801VX3    +(  -4.92691VX5 

+(  -0.13851VX8    +(  -0.13851VX9    +(  0.00061FC81 

+(  0.00261FCB3+I  22 . 4889 1Y1 1    + (  -75.21281Y13 

+(  22.48891Y21    +(  -75 .2128 1 Y2 3   + (  0.12371Y011 

+  (  -0.61901YC13+(  0.12371YD21+I  -0.6190IYD23 


ECUATICN 


1.0    AX6    =+( 

2.93181X2      +( 

384.4605  1X4       +( 

-42736.91021X6 

+  ( 

69.04131X7       +( 

-31.C9571X8      +1 

31.09571X9 

+  ( 

42.40441X10    +( 

-42.40441X11    +( 

5389.75001X14 

+  1 

-5389.75001X15    +( 

-6699.26521X22    +( 

-6699.28521X23 

+  ( 

20.41931X26   +( 

20.41931X27    +( 

-32.28801X28 

+  ( 

-32.28801X29    +{ 

0.02961VX2    +( 

0.10691VX4 

+  ( 

-1.68451VX6    +( 

0.16131VX7    +( 

-0.17101VX3 

+  ( 

0.1710IVX9    +( 

-0.0868 1VX26+1 

-0.08O81VX27 

♦  ( 

-0.C8681VX28+ I 

-0.08681VX29+( 

-0.0015 1FDB2 

+  1 

27.76401Y11    +( 

30.82201Y12    +( 

-27.76401Y21 

+  ( 

30.82201Y22    +( 

-16.14401Y32    +( 

-16.14401Y42 

+  ( 

-16.14401Y72    +( 

-16.14401Y82    +( 

0.15271YC11 

+  1 

0.07201YC12+( 

-0.15271YD21+C 

0.07201YC22 

+  ( 

-0.04341YC3  2+( 

-0.04341YD42+I 

-0.04341YC72 

+  ( 

-0.04341Y082 

ECUATION           7 

1.0    AX7   =+( 

742.30001X2      +( 

-1150.56421X4       +( 

2724.24001X6 

+  i 

-9302.76911X7      +( 

1.73431VX2    +( 

-2.68821VX4 

+  ( 

6.36501VX6    +( 

-71.66301VX7    +( 

-371.1499IY12 

+  ( 

-1721.8335IY13   +( 

-371.14991Y22    +( 

1721.83351Y23 

+  ( 

-0.8672)YC12+( 

-14.1708)Y013+( 

-0.36721YD22 

+  ( 

14.17081YD23 

ECUATICN           8 

1.0    AX8    =+( 

-15719.30081X1       +( 

-24364.90231X5       +( 

-38669.46881X6 

+  ( 

-17605. 6C941X8      +( 

-86.45621VX1    +{ 

-134.00701VX5 

+  ( 

-212.6821IVX6    +1 

-96.d309)VX8    +( 

15719. 30081Y11 

+  ( 

86.45621Y011 
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E0UATION         18 
1.0    AX18»+(         23389.24221X5      +(         -1588.70311X14   +  <       -57752.27341X18 
+(  31185.65231X20    +(  -1 1. 8270 1 VX14+ <  -23. 6540 1 VX1 3 

+(  1588.70311Y31    +(  11.82701Y031 

EQUATION         19 
1.0   AX19-H  23389.24221X5      +(         -1586. 7C31 1 X15   +  <       -57  7  52.27341X19 

+(         31185.65231X21    +(  -11. 8270 1 VX 15* <  -23. 65401 VX1 S 

+  (  1588.70311Y41    <■(  ll.S270)YD41 

EQUATION         20 
1.0   AX20=-H         -1588.7C311X16    +(         31185.6523  1X18    *■(       -34363.05861X20 
+(  -11.8270JVX16+1  -23.6540JVX20+1  1588. 70' 1) Y71 

+<  11.8270IYD71 


EOJATICN         21 
1.0   AX21  =  +  (         -1588.70311X17    ♦  I         3118  5.65  231X19    +(      -34363.05861X21 
♦(  -11.82701VX17+1  -23.6540IVX21+I  1568. 7031) Y81 

+(  11.82701YD81 


EQUATION         22 
1.0  AX22=-H    -234187.87501X1       * (    -1 56905. 8 125 1 X5      +(-838245.68751X6 
+  <    -216580. 25C01X14    +(       450768. 25C0 1 X16    +(    -935  251.18751X22 
+(    -408571. OC001X24 


ECUATICN         23 
1.0   AX23=+<       234187.87501X1       +(       156905.81251X5       ♦<    -838245.66751X6 
+  (       216580. 25C01X15    +<    -450768. 25  CO  I  X 17   +  <    -935251.18751X23 
+(    -406571. OC001X25 


EQUATION  24 

1.0   AX24-+I    -450767.81251X14    +(       450767.81251X16    +(    -4085  71.00001X22 
+(    -565087.37501X24 


ECUATICN         25 
1.0   AX25=  +  <       450767.81251X15    *•(    -450  767.81251X17    +(    -408571.00001X23 
+(    -565087.37501X25 


,     „  ECUATION          26 

1.0    AX26=  +  (  210C. 6174-1X2       ♦<       2341S7.6  7  501X3       +(    -839653.06251X4 

+(  545657.81251X5       +1          -4694.43751X6       + (       216  580.25001X10 

+(  -450768.25001X12    -M    -939452.43  751X26    +(    -4065 71 . OCOO) X2d 

+<  5.64671VX2    +(                  -3.78461VX4    +(               -13.1615IVX6 

+(  -11.2974)VX26t(         -2100. 6 1 74 ) Y32    +(                 -5.64871YD32 


,  „  .„,,  EQUATION    27 

1.0  AX27=+(     2100.61741X2   <■  (  -234167.87  501X3   +  (  -839653.06251X4 


Ill 


♦(  -545657.81251X5      +(         -4894.43751X6  ♦(  -216  580.25001X11 

+  (  450768.25001X13    +(    -939452.43751X27  * I  -4085 7  I. OC001 X29 

+  (  5.64871VX2    +(                 -3.76461VX4  +(              -13.16151VX6 

+  (  -U.29741VX27+I         -2100. 6174) Y42  +  1                 -5.64871Y042 


EOUATICN        28 

1.0    AX28=+<           2100.61741X2      * (         -1407.4133  1X4      +(  -4894.43  751X6 

+  (       450767.81251X10    *•  (    -450767.8125  1X12    +1  -408  570.62501X^6 

+(    -569288.62501X28    +(                     5.64671VX2-H  -3. 7846) VX4 

+  {               -13.16151VX6    ♦(               -U.2974)VX28+(  -2100.  61 74)  Y72 
+(                  -5.64871Y072 

ECUATION  29 

1.0    AX29=+<             2100.6174JX2      +(          -1407.41331X4      +  <  -4894.43  751X6 

+  (    -450767.81251XU    *  (       450767.81251X13   +(  -408  570.62501X27 

+  (    -569286.62501X29+1                    5.6487)VX2+(  -3.7846IVX4 

<■(               -13.16151VX6   <■(               -11.2974JVX29+1  -2 100.  6174)  Y82 
+(                -5.64871YD82 
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A  dual  rear  wheeled  farm  tractor  is  modeled  as  eight  lumped  masses 
with  elastic  rear  axles.   The  equations  of  motion  are  derived  using  Newton's 
laws,  and  elementary  beam  theory.  These  equations  are  checked  by  deriving 
the  same  equations  using  the  energy  method  with  the  Lagrange  equation  for 
non-conservative  systems.  The  equations  of  motion  are  integrated  numerically 
using  the  Runge-Kutta-Gill  method,  to  give  the  motion  of  the  vehicle  as  a 
function  of  time. 

Computer  programs  are  presented  which  generate  the  equations  of  motion, 
using  either  Newton's  laws  or  the  energy  method,  for  a  given  set  of  tractor 
parameters.  A  program  to  perform  the  numerical  integration  is  also  pre- 
sented. This  program  includes  provisions,  by  which  the  system  may  be  forced. 
The  forcing  may  be  due  to  a  drawbar  load  or  surface  profile. 

The  computer  programs  are  demonstrated  by  simulating  the  motion  of  a 
tractor  traversing  a  half  sinusoid  bump  with  the  left  rear  wheels.  The 
results  of  a  system  with  soft  rear  axles  are  compared  with  the  results  of 
a  system  with  stiff  axles. 


